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, Abstract 

With its rapid-response capability and multiwavelength complement of instruments, the Swift 

satellite has transformed our physical understanding of 7-ray bursts (GRBs). Providing high- 

^ quality observations of hundreds of bursts, and facilitating a wide range of follow-up observations 

within seconds of each event. Swift has revealed an unforeseen richness in observed burst prop- 

erties, shed light on the nature of short-duration bursts, and helped realize the promise of GRBs 

^— ( as probes of the processes and environments of star formation out to the earliest cosmic epochs. 

These advances have opened new perspectives on the nature and properties of burst central 

engines, interactions with the burst environment from microparsec to gigaparsec scales, and the 

possibilities for nonphotonic signatures. Our understanding of these extreme cosmic sources has 

thus advanced substantially; yet, more than 40 years after their discovery, GRBs continue to 

^ present major challenges on both observational and theoretical fronts. 
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1 INTRODUCTION 
1.1 Setting the Stage 

Gamma-ray bursts (GRBs) are among the most fascinating phenomena in the 
Universe. They are bright flashes of radiation with spectral energy distributions 
peaking in the 7-ray band. They have durations measured in seconds and appear 
to be capable of producing directed flows of relativistic matter with kinetic lumi- 
nosities exceeding 10^^ erg s~^, making them the most luminous events known. 
All evidence points to a gravitational power source associated with the cata- 
clysmic formation of a relativistic star or to a precursor stage whose inevitable 
end point is a stellar mass black hole. 

The field of GRB astronomy has been greatly stimulated by the launch of the 

in 2004 with its rapid response and panchro- 
matic suite of instruments onboard, and by the development of new-technology 
robotic telescopes on the ground. A multidisciplinary approach is now emerg- 
ing with data combined across the electromagnetic spectrum to learn about the 
physical processes at play; "spectral chauvinism" can no longer be tolerated in 
the modern study of GRBs. Even non-photonic neutrino and gravitational wave 
instruments are becoming more sensitive and may soon be detecting signatures 
related to GRBs. 

While interesting on their own, GRBs are now rapidly becoming powerful tools 
to study detailed properties of the galaxies in which they are embedded and of the 
universe in general. Their apparent association with massive star formation and 
their brilliant luminosities make them unique probes of the high-redshift universe 
and galaxy evolution. Absorption spectroscopy of GRB afterglows is being used 
to study the ISM in evolving galaxies, complementary to the traditional studies 



Swift sateinte (Gehrels et aL]|2004 



2 



3 



of quasar absorption line systems. Possibly the most interesting use of GRBs in 
cosmology is as probes of the early phases of star and galaxy formation, and the 
resulting reionization of the universe at z ~ 6 — 20. GRBs are bright enough 
to be detectable, in principle, out to much larger distances than those of the 
most luminous quasars or galaxies detected at present. Thus, promptly localized 
GRBs could serve as beacons which, shining through the pregalactic gas, provide 
information about much earlier epochs in the history of the universe. 



Before the advent of Swift, the study of GRBs had evolved somewhat unsys- 
tematically, and, as a result, the field has a great deal of historical curiosities 
such as complex classifications schemes which are now becoming streamlined as 
the field matures. Objects which were once thought to be different are now 
found to be related and the style of research has shifted from piecewise studies 
to a more general statistical approach. Although leaps in understanding can still 
come from extraordinary events as we show in several examples in this review, 
the applications to broader astrophysics are coming from the compilations of 
hundreds of events. The literature on this subject has therefore become quite 
large, and we apologize for referring now and then only to the most recent com- 
prehensive article in a given topic. There are several recent summary articles 
that give excellent reviews in specific areas related to GRBs. These include the 



supernova-burst connection 


(Woosley and Bloom 


2006 


1, short GRBs ( 


Lee and 


Ramirez-Ruiz 2007 


Nakar 2007a|), afterglows (van Paradijs et al. 2000 


Zhang 


2007 


) and theory ( 


Meszaros 


2002). Our objective 


bere is to summarize the field 



of GRB astronomy, from the Swift era and prior to the next steps with the Fermi 
Gamma Ray Observatory (Atwood et al. 2009), interpreting past findings while 
looking ahead to future capabilities and potential breakthroughs. 



1.2 A Burst of Progress 

The first sighting of a GRB came on July 2, 1967, from the military Vela satellites 
monitoring for nuclear explosions in verification of the Nuclear Test Ban Treaty 



(Klebesadel et al. 1973). These 7-ray flashes, fortunately, proved to be different 
from the man-made explosions that the satellites were designed to detect and a 
new field of astrophysics was born. Over the next 30 years, hundreds of GRB 
detections were made. Frustratingly, they continued to vanish too soon to get 
an accurate angular position for follow-up observations. The reason for this is 
that 7-rays are notoriously hard to focus, so 7-ray images are generally not very 
sharp. 

Before 1997, most of what we knew about GRBs was based on observations 
from the Burst and Transient Source Experiment (BATSE) on board the Comp- 



ton Gamma-Ray Observatory, whose results have been summarized in Preece 



et al. (2000). BATSE, which measured about 3000 events, revealed that between 
two and three visible bursts occur somewhere in the universe on a typical day. 
While they are on, they can outshine every other source in the 7-ray sky, in- 
cluding the Sun. Although each is unique, the bursts fall into one of two rough 
categories. Bursts that last less than two seconds are short, and those that last 
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longer - the majority are long. The two categories differ spectroscopically, with 
short bursts having relatively more high-energy 7-rays than long bursts do. 

Arguably the most important result from BATSE concerned the distribution 
of bursts. They occur isotropically - that is, evenly over the entire sky — sug- 
gesting a cosmological distribution with no dipole and quadrupole components. 
This finding cast doubt on the prevailing wisdom, which held that bursts came 
from sources within the Milky Way. Unfortunately, 7-rays alone did not provide 
enough information to settle the question for sure. The detection of radiation 
from bursts at other wavelengths would turn out to be essential. Visible light, 
for example, could reveal the galaxies in which the bursts took place, allowing 
their distances to be measured. Attempts were made to detect these burst coun- 
terparts, but they proved fruitless. 

A watershed event occurred in 1997, when the BeppoSAX satellite succeeded 
in obtaining high-resolution X-ray images (Piro et al. 1999) of the predicted 
fading afterglow of GRB 970228 - so named because it occurred on February 28, 
1997. This detection, followed by a number of others at an approximate rate of 
10 per year, led to positions accurate to about an arc minute, which allowed the 
detection and follow-up of the afterglows at optical and longer wavelength^ (e.g.. 



van Paradijs et al. 1997). This paved the way for the measurement of redshift 



distances, the identification of candidate host galaxies, and the confirmation that 
they were at cosmological distances (Metzger et al. 1997). 



Among the first GRBs pinpointed by BeppoSAXwas GRB 970508 (Metzger 



et al. 1997). Radio observations of its afterglow provided an essential clue. The 



glow varied erratically by roughly a factor of two during the first three weeks, 
after which it stabilised and then began to diminish (Frail et al. 1997). The 
large variations probably had nothing to do with the burst source itself; rather 
they involved the propagation of the afterglow light through space. Just as the 
Earth's atmosphere causes visible starlight to twinkle, interstellar plasma causes 
radio waves to scintillate. Therefore, if GRB 970508 was scintillating at radio 
wavelengths and then stopped, its source must have grown from a mere point to 
a discernible disk. "Discernible" here means a few light-weeks across. To reach 
this size the source must have been expanding at a considerable rate ~ close to 
the speed of light (Waxman et al. 1998). 

The observational basis for a connection between GRBs and supernovae was 
prompted by the discovery that the BeppoSAX error box of GRB 980425 con- 
tained supernova SN 1998bw (Galama et al. 1998). A number of other GRBs 
have since shown a 1998bw-like temporal component superposed on the power 
law optical lightcurve (Woosley and Bloom 2006), but they still lacked a clear 
spectroscopic detection of an underlying supernova. Detection of such signature 



came with the discovery by the HETE-2 mission of GRB 030329 (Stanek et al. 



2003 , Hjorth et al.||200"3 ) . Due to its extreme brightness and slow decay, spectro- 



scopic observations were extensive. The early spectra consisted of a power-law 
decay continuum {F^ oc typical of GRB afterglows with narrow emission 



^We note, however, that the first optical afterglow detection of GRB 970228 (van Paradijs 
et al.|1997 ) was based on the X-ray prompt detection by BeppoSAX. 
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features identifiable as Ha, [OIII], H/3 and [Oil] at z = 0.1687 (Kawabata et al. 



2003 Matheson et al.|[2003[ ) making GRB 030329 the second nearest burst over- 
all at the time and the nearest classical burslj^ A major contribution to our 
understanding of X-ray prompt emission came also from the HETE-2 mission 
(Lamb et al. 2004 1 which was active from 2000 to 2006. Dozens of bursts in the 
"X-ray flash" category were observed and were found to be similar in origin to 
the classical long GRBs (Matsuoka et al. 2004). 

Swift is the current GRB discovery mission. It is a space robot designed specif- 
ically with GRBs in mind. It combines a wide-field hard X-ray burst detection 
telescope (Burst Alert Telescope - BAT; Barthelmy et al. 2005a I, with narrow- 
field X-ray (X-Ray Telescope - XRT; [Burrows et al.|[2005a| , and UV-optical (UV 
Optical Telescope - UVOT; Roming et al. 2005) telescopes. A powerful and fast 
on-board burst detection algorithm ( Fenimore et al.||2003 | provides the burst co- 
ordinates to the spacecraft, which autonomously repoints the observatory so that 
X-ray and optical observations typically commence within two minutes of the 
burst trigger. The mission was designed to find counterparts for all burst types, 
including the previously-elusive short GRBs. Burst positions and other data are 
provided promptly to the ground for ground observers. 

The burst detection rate for Swift is about 100 GRBs per year, resulting in a 
current data set as of December 2008 of 380 bursts. Of these, there are 126 with 
redshift determination, mostly from spectrographs on large optical telescopes and 
new robotic telescopes on the ground. These now far outnumber the ~ 40 GRB 
redshifts available prior to Swift. More than 95% of the Swift bursts have X-ray 
afterglow detection and ~ 60% have optical afterglows (UVOT -|- ground). To 
date 33 short bursts have been localized with 8 having redshift determinations. 
The new data have enabled much more detailed studies of the burst environment, 
the host galaxy, and the intergalactic medium between galaxies. Swift and follow- 
up observations have also transformed our view of GRB sources. For example, as 
discussed in Section [3) the old concept of a sudden release of energy concentrated 
in a few brief seconds has been discarded. Indeed, even the term "afterglow" 
is now recognized as misleading - the energy radiated during both phases is 
comparable. 

Our primary intention in this review is to describe the most important ob- 
servational discoveries in the Swift era as well as to explain some of the new 
understanding of key mechanisms that are believed to operate in these objects to 
an astronomical audience with little prior exposure to GRBs. Four sections fol- 
low. Section |2] provides a description of our current knowledge of what constitutes 
a GRB. Section [3] is a summary of the observations of the prompt and afterglow 
emission and Section IV of the observations of host galaxies and progenitor clues. 
In Section V, we examine our current progress in our understanding of the basic 
physical processes at work. Section VI is a look forward at future prospect for 
GRB study. 



The other GRBs with supernova associations have been underluminous events (Kaneko et al. 



20071. 
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2 WHAT IS A GAMMA-RAY BURST? 



GRBs are sudden, intense flashes of 7— rays which, for a few Winding seconds, 
Hght up an otherwise fairly faint 7— ray sky. Spectra extending over many decades 
in photon energy have now been measured for hundreds of GRBs. In Figure [l] 
representative spectra are plotted in the conventional coordinates v and vF^^ the 
energy radiated per logarithmic (natural log) frequency interval. Some obvious 
points should be emphasized. First, we measure directly only the energy radiated 
in the direction of the earth per second per steradian per logarithmic frequency 
interval by a source. The apparent bolometric luminosity may be quite different 
from the true bolometric luminosity if the source is not emitting isotropically. 
Second, there is striking evidence for a characteristic photon energy (peak in 
the vFy spectrum), which appears to be related to the overall spectral luminos- 
ity normalization. In contrast, the spectra of many galactic and extragalactic 
accretion systems are often well fltted by single power-laws. A simple power- 
law contains little information, whereas a complex spectrum composed of many 
broken power-laws tell us much more, as each break frequency must be explained. 

At cosmological distances the observed GRB fluxes imply energies that can 
exceed 10^^(O/47r) erg, where is the solid angle of the emitting region (Figure[2] 
see also [Bloom et ah 20011 . This is the mass equivalent of O.OGMq for the 
isotropic case. Compared with the size of the sun, the seat of this activity is 
extraordinarily compact, with sizes of less than milli- light-seconds (<300 km) 



as indicated by rapid variability of the radiation flux (Bhat et al. 1992 1 . It is 



unlikely that mass can be converted into energy with better than a few (up to 
ten) percent efficiency; therefore, the more powerful GRB sources must "process" 
upwards of lO^^(r2/47r)-/Vf0 through a region which is not much larger than the 
size of a neutron star (NS) or a stellar mass black hole (BH). No other entity 
can convert mass to energy with such a high efficiency, or within such a small 
volume. 

The observed 7-rays have a nonthermal spectrum. Moreover, they commonly 
extend to energies above 1 MeV, the pair production threshold in the rest frame. 
These facts together imply that the emitting region must be relativistically ex- 
panding CGuilbert et al. 1983, Goodman 1986 Paczyhski 19861. We draw this 



conclusion for two reasons. First, if the region were indeed only a light second 
across or less, as would be implied by the observed rapid variability in the ab- 
sence of relativistic effects, the total mass of baryons in the region would need to 
be below about 1O~^^M0 in order that the electrons associated with the baryons 



should not provide a large opacity (Piran and Shemi 1993 Paczyhski 1990). 



Second, larger source dimensions are required in order to avoid opacity due to 
photon-photon collisions. If the emitting region is expanding relativistically, then, 
for a given observed variation timescale, the dimension R can be increased by F^. 
The opacity to electrons and pairs is then reduced by F^, and the threshold for 
pair production, in the observer frame, goes up by ~ F from its rest frame value 
(Fenimore et al. 1993, Woods and Loeb 1995 Baring and Harding||1997 Granot 



et al. 20081. Best-guess numbers are Lorentz factors F in the range 10 to 10 
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energy [keV] 



Figure 1 Gamma-rays are excellent probes of the most energetic phenomena in na- 
ture, that typically involve dynamical non-thermal processes and include interactions 
of high energy electrons with matter, photons and magnetic fields; high energy nuclear 
interactions; matter- antimatter annihilation and possibly other fundamental particle in- 



teractions. Here shown are representative spectra vF,^ cx v^N^v) of GRBs (Kaneko et al. 
2007 Kaneko et al.|2008 ) together with the Crab pulsar nebula ( Kuiper et al.|2001 1 and 
the galactic black hole candidate Cygnus X-1 ( [McConneU eraI1[2002l ). 



(Lithwick and Sari 2001), allowing rapidly- variable emission to occur at radii in 
the range 10^^ to 10^"^ cm. 

Because the emitting region must be several powers of ten larger than the 
compact object that acts as a trigger, there are further physical requirements. 
The original internal energy contained in the radiation and pairs would, after ex- 
pansion, be transformed into relativistic kinetic energy. A variant that has also 
been suggested is based on the possibility that a fraction of the energy is carried 



by Poynting flux (Blandford and Znajek 1977 Usov 1992). This energy cannot 



be efficiently radiated as 7-rays unless it is re- randomized (Meszaros et al.||1994 


Narayan et al. 


1992 


Rees and Meszaros 


1994 


Paczyhski and Xu 


1994 


). Impact 


on an external medium (or an intense external radiation field; e. 


g., Shaviv and 



Dar 1995 ) would randomize half of the initial energy merely by reducing the ex- 
pansion Lorentz factor by a factor of 2. For an approximately smooth distribution 
of external matter, the bulk Lorentz factor of the fireball thereafter decreases as 
an inverse power of the time. In the presence of turbulent magnetic fields built 
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Figure 2 Apparent isotropic 7-ray energy as a function of redshift and observed duration. 
The energy is calculated assuming isotropic emission in a common comoving bandpass 
for a sample of short and long GRBs with measured redshifts. This spread in the inferred 
luminosities obtained under the assumption of isotropic emission may be reduced if most 
GRB outflows are jet-like. A beamed jet would alleviate the energy requirements, and 
some observational evidence does suggest the presence of a jet. 



up behind the shocks (Rees and Meszaros 19921, the electrons produce a syn- 



chrotron power-law radiation spectrum which softens in time, as the synchrotron 
peak corresponding to the minimum Lorentz factor and field decreases during 



the deceleration (Katz 1994 Sari et al. 1996). Thus, the GRB radiation, which 



started out concentrated in the 7-ray range during the burst, is expected to pro- 
gressively evolve into an afterglow radiation which peaks in the X-rays, then UV, 
optical, IR and radio (Figure [3|. Detailed predictions (Meszaros and Rees 19971 
of the afterglow properties, made in advance of the observations, agreed well 
with subsequent detections at these photon energies, followed up over periods 
of up to months. The detection of diffr active scintillation in the radio afterglow 
of GRB 970508, provided the first determination of the source size and a direct 
confirmation of relativistic source expansion dFrail et al.|T997l ), which was further 
strengthened by the size measurement of the afterglow image of GRB 030329 by 
radio interferometry with the VLBA ( Taylor et al.||2004 ). 

The complex time-structure of some bursts suggests that the central engine 
may remain active for up to 100 s ( Ramirez-Ruiz and Fenimore||2000 1 or possibly 
longer (Falcone et al. 2007). However, at much later times all memory of the 
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Figure 3 The evolving synchrotron afterglow of a 7— ray burst. Shown is a theoretical 
model ( |Gou et al.|[2007[ ) for the afterglow of the Swift GRB 050904. The model is pre- 
sented without extinction and as it would have been observed at redshift z = 2; the burst 
itself occurred at z = 6.29. The evolution of the synchrotron peak to lower frequencies 
is clearly visible. More subtle effects including evolution of the synchrotron cooling and 
self-absorption frequencies, and the associated synchrotron self-Compton emission of the 
blastwave at higher frequencies, are not readily visible in this model. 



initial time-structure would be lost: essentially all that matters is how much 
energy and momentum has been injected, its distribution in angle and velocity. 
We can at present only infer the energy per solid angle, but there are reasons 
to suspect that bursts are far from isotropic. Due to relativistic beaming, an 
observer will receive most emission from those portions of a GRB blast wave that 
are within an angle ~ l/F of the direction to the observer. The afterglow is thus 
a signature of the geometry of the ejecta - at late stages, if the outflow is beamed, 
we expect a spherically-symmetric assumption to be inadequate; the deviations 
from the predictions of such a model would then tell us about the ejection in 
directions away from our line of sight (Rhoads 19991. 

The appearance of achromatic breaks in the development of GRB afterglows 
has been interpreted as indicating that they are jet flows beamed towards us. 
Collimation factors of Qi/47r ^ 0.01 (corresponding to half opening angles of ^ 8 



degrees) have been derived from such steepening (Frail et al. 2001 Bloom et al. 



2003). If GRB sources are beamed, then this reduces the energy per burst by two 



or three orders of magnitude at the expense of increasing their overall frequency. 

As regards the central engine trigger, there remain a number of key questions. 
What are the progenitors? What is the nature of the triggering mechanism, the 
transport of the energy and the time scales involved? Does it involve a hyper- 
accreting compact object? If so, can we decide between the various alternative 
ways of forming it? The presence of some GRBs (in the short duration category) 
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in old stellar populations rules out a source uniquely associated with recent star 



formation and, in particular, massive star origin for all bursts (Gehrels et al. 



2005 , Bloom et al.||2006 ). An understanding of the nature of these sources is thus 



inextricably linked to the metabolic pathways through which gravity, spin, and 
energy can combine to form collimated, ultrarelativistic outflows. These threads 
are few and fragile, and the tapestry is as yet a poor image of the real universe. 
If we are to improve our picture-making we must make more and stronger ties to 
physical theory. But in reconstructing the creature, we must be guided by our 
eyes and their extensions. The following sections provide a detail summary of 
the observed properties of these ultra-energetic phenomena. These threads will 
be woven in Section |5l 



3 BURST AND AFTERGLOW OBSERVATIONS 

The most direct diagnostics of the conditions within GRBs come from the ra- 
diation we observe, which we summarize in this section. There are many ways 
to organize a discussion of the properties of the radiation emanating from GRB 
sources. We do not intend in this article to a give a detailed review of individual 
events since there are now sufficiently many examples that we are likely to be 
led seriously astray if we test our theories against individual events. For this 
reason, we center our discussion on major trends even in cases in which the gen- 
eralizations we describe are based on data that do not yet have the statistical 
weight of a "complete" sample. It should be noted that there are inherent biases 
in the discovery of a GRB at a given redshift which are often difficult to quan- 
tify, such as complex trigger efficiencies and non-detections. Continued advances 
in the observations will surely yield unexpected revisions and additions in our 
understanding of the properties of GRBs. 



3.1 Prompt High-Energy Emission 

3.1.1 Taxonomy. The manifestations of GRB activity are extremely di- 
verse. GRBs are observed throughout the electromagnetic spectrum, from GHz 
radio waves to GeV 7-rays, but until recently, they were known predominantly 
as bursts of 7-rays, largely devoid of any observable traces at any other wave- 
lengths. 7-ray properties provide only one of several criteria for classifying GRB 
sources. Part of the problem is observational because it is not possible to obtain 
full spectral coverage in all objects and it is not easy to reconcile a classification 
based on host galaxy properties with one based on the prompt 7-ray properties. 
The major impediment to serious taxonomy is more fundamental. GRBs are het- 
erogeneous objects, especially in their directly observed properties. The success 
of a classification scheme, we believe, should be measured by the extent to which 
newly recognized properties distinguish subsets defined by differences in other 
properties. By this criterion, the taxonomy of GRBs has met with only mixed 
success. As new non-7-ray selection techniques are introduced (e.g. age of stellar 
populations in host galaxies or the presence of type Ic supernova signatures), the 
class boundaries (e.g. short and long duration events) have blurred where the 
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defined subclasses transcend traditional boundaries. On the other hand, many 
new properties do correlate with old ones. This is all the more remarkable in 
that the conventional diagnostics (e.g. burst duration) measure properties on 
scales several orders of magnitude larger than that which we believe to be the 
characteristic of the powerhouse. 



3.1.2 Observed durations and redshifts. GRB traditionally have been 
assigned to different classes based on their duration - usually defined by the time 
during which the middle 50% (T50) or 90% (T90) of the counts above background 
are measured. On the basis of this criterion, there are two classes of GRBs - 
short and long - where ~2 s duration separates them. The initial hints for the 
existence of such classes ( Cline and Desai|1974 Mazets et al.|198T |, were followed 



by stronger evidence from ISEE-3 and Konus-WIND data (Norris et al. 



and definite proof using large statistics from BATSE (Kouveliotou et al. 



19841 



19931. 



BATSE results also showed that short bursts have a harder spectrum than long 
bursts (Kouveliotou et al. 19931 although this is less prominent in observations 
by Konus-WIND, HETE-2, and Swift ( [Sakamoto e t al.||2006[). 

The duration and redshift distributions for Swift GRBs are shown in Figure |4j 
The blue histogram in the left panel is for observed durations while the orange 
histogram shows the durations corrected to the source frame Tgo/{l + z) for those 
bursts with redshift determinations. In the source frame, the typical long burst 
duration is ~ 20s compared to ~ 50s in the observer frame. Swift has been 
detecting a lower fraction (~10%) of short bursts than BATSE did (25%). This 
is because Swift observes in a softer energy band (15 — 150 keV) than BATSE 
(50 keV — 2 MeV) and because Swift requires a sky image of the event for 
burst detection and the image part of the trigger algorithm is less sensitive to 
short bursts due to their lower fluences. Shown in Figure |4] (ri^/it panel) are the 
measured redshift distributions. The blue histogram is for Swift events while the 
gray one is for pve-Swift bursts. As clearly seen. Swift is currently detecting GRBs 
at a higher average redshift: < z >~ 2.5 for Swift bursts while < z >~ 1.2 for 
pie-Swift events. The reason for this difference is the higher sensitivity of Swift 
compared to BeppoSAX and HETE-2. 

GRBs have also been classified according to their spectral properties, albeit 
less successfully. In particular, bursts with lower spectral energy peaks (E'peak) 
have been denoted X-Ray Flashes (XRFs) based on observations by BeppoSAX, 
BATSE and HETE-2 ( [Heise et al.]|200T| [Barraud et al][2003l [Kippen et aT][2003 



Sakamoto et al.||2005 ) . These events are closely related to common long duration 



GRBs and appear to form a continuum of all parameters between the two types 
with no striking evidence for a discerning characteristic (Granot et al. 2005b I . 



3.1.3 Observed correlations. There is a great deal of diversity in the 7-ray 
prompt lightcurves of GRBs. Both long and short bursts can have temporal pro- 
files ranging from smooth, single-peaked pulses to highly structured multi-pulses. 
The prompt emission can be characterized by a variety of spectral and temporal 
parameters which include duration, variability, lag, pulse rise/fall time, fluence. 
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Figure 4 Duration and redshift distributions for Swift GRBs. Left panel shows the 
duration distribution. The blue histogram is the measured Tgo distribution while the 
orange one is corrected to the source frame: 190/(1 + 2). Right panel shows the redshift 
distribution for Swift GRBs in blue and pre- Swift GRBs in gray. Swift is detecting higher 
redshift bursts on average than pie-Swift. The solid (broad) theory curve illustrates the 
evolution of a comoving volume element of the universe while the dotted (narrow) curve 
is a convolution of the comoving volume with a model for the star formation rate as 



calculated by Porciani and Madau ( 2001 1 



£■180 and Epeak. A schematic diagram illustrating the most widely discussed 7- 
ray prompt correlations is shown in Figure |5] with detailed references given in its 
caption. These correlations are often based upon statistical analysis of quantities 
whose physical causes are poorly understood, but almost certainly depend on 
many variables. Interpretations of these correlations must therefore be done with 
caution. 

The prompt GRB lightcurves can generally be dissembled into a superposition 



of individual pulses as described by Norris et al. (19961 with rise times shorter on 



average than decay times (panel f). The variability or spikeness of the lightcurve 
is found to be correlated with peak luminosity or total isotropic energy of the 
burst (panel a). The time lag of individual peaks seen at different energy bands 
is observed to be anticorrelated with luminosity for long bursts (panel b). For 
short bursts, the lag is small or not measurable. The Ep^ak is also found to be 
correlated with Eiso for long bursts including XRFs with short bursts as clear 
outliers {panel c) . The total isotropic energy emission is correlated with duration 
{panel d), with short and long bursts on approximately the same correlation line, 
albeit with a wide spread. Short bursts detected by Swift have lower E'iso on 
average than long bursts. There is a group of outliers belonging to the long burst 
category which are characterized by being significantly underluminous. These 
are GRBs: 980425, 031203, 060614. 

Numerous researchers have studied ways to determine the absolute luminosity 
of a GRB using correlations such as those illustrated in Figure |5j These include 
the lag, variability and -Bpeak correlations discussed above. Other interesting cor- 
relations have included -Epeak versus E^ (emitted energy corrected for beaming. 



Ghirlanda et al. 2004) and -Epeak versus a duration corrected peak luminosity 
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(Firmani et al. 2006). The goal is to derive a method to determine the burst 



luminosity independently of a redshift distance determination, thus attempting 
to make GRBs standard candles which could be used, in principle, to deter- 
mine the cosmic expansion history of the universe to higher redshift than with 



supernovae. Although such efforts are currently underway (Schaefer 2006, and 



references therein) , it is not clear at present whether any of these correlations are 
tight enough for significant progress to be made ( Bloom et al.||200^ |. 



3.1.4 Soft 7 repeaters and short bursts. It has been noted (Hurley 



et al. 2005 Palmer et al. 2005 Nakar et al. 2006b I that the giant flare (GF) 



observed from the putative galactic magnetar source SGR1806-20 in December 



2004 (Gaensler et al. 20051 could have looked like a typical short GRB had it 



occurred much farther away, thus making the tell-tale periodic signal character- 
istic of the neutron star rotation in the fading emission undetectable. The two 
previously recorded GFs of this type, one each from SGR 0520—66 on 5 March 
1979 dFenimore et"al]|1996a[ ) and SGR 1900-^14 on 27 August 1998 ( [Hurley et alT 
19991, would have been detectable by existing instruments only out to ~ 8 Mpc, 



and it was therefore not previously thought that they could be the source of short 
GRBs. The main spike of the 27 December event would have resembled a short, 
hard GRB if it had occurred within ~ 40 Mpc, a distance scale encompassing the 
Virgo cluster ( Palmer et al.|[2005 |. However, the paucity of observed GFs in our 
own Galaxy has so far precluded observationally based determinations of either 
their luminosity function or their rate. The observed isotropic distribution of 
short BATSE GRBs on the sky and the lack of excess events from the direction 
of the Virgo cluster suggests that only a small fraction, < 5%, of these events 
can be SGR GFs within 40 Mpc ( [Palmer et al]|2005[ ). 

Before Swift detected short GRBs and their associated afterglow signatures, 
searches for nearby galaxies within narrow Inter Planetary network (IPN) error 
boxes revealed that only up to ~ 15% of them could be accounted for by SGRs 
capable of producing GFs ( Nakar et al.|[2006b ). A recent, intriguing candidate is 
short GRB 070201 which was observed by the IPN to have a location consistent 



with the arms of the nearby (0.8 Mpc) M31 galaxy ( [Mazets et al.[[2008[ ). A 
LIGO search for gravitational waves ( [Abbott et al.|2008 | at the time of the burst 
turned up with no signal, thus excluding a compact merger origin. If the GRB 
was really in M31, it may have been an SGR GF. While the fraction of SGR 
events among what are now classified as short GRBs may not be dominant, it 
should be detectable and can be tested with future Swift observations. It is also 



worth noting that some short GRBs likely originate in the local univserse (Tanvir 
eral][2005|). 



3.2 Afterglow Observations 

3.2.1 X-ray observations. Swift was designed to investigate the GRB af- 
terglows by filling the temporal gap between observations of the prompt emission 
and the afterglow. The combined power of the BAT and XRT has revealed 
that prompt X-ray emission smoothly transitions into the decaying afterglow 
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Figure 5 Scliematic diagrams illustrating the most widely discussed correlations between 
various prompt emission properties for long (L), short (S ), and underluminous (UL) 
GRBs. (a) Variability scaled to the burst frame versus Ej^n (iFenimore and Ramirez-Ruiz 



2000 Reichart et al. 



ot the lightcurve and 



20011 |SchaefeF||2006[ ). The variability is a measure of the spikiness 
I IS dehned as the mean square of the time signal after removing 
(b) Spectral lag scaled to the burst frame versus peak 
2006J. (c) -Epcak scaled to the burst 



low frequen cies by s moothing, 

luminosity (Norris and BonncU 2006: Gchrcls ct al ^ 

frame versus Ej^p (|Amati et a l. 2002 for BeppoSAX GK^s; |AmalTj|2006| for Swift GRBs; 
Lloyd- Ronning and Kamirez-Ruiz,|2002| for BATSE events), (d) Duration scaled to the 
burst fram e versus Ejso. (e) Spectral hardness versu s observed duratio n (|Kouveliotou 
eFallpMI . (f) Pulse rise time versus its decay time ( [Norris et al.||1996| . 




Figure 6 Representative examples of X-ray afterglows of long and short Swift events with 
steep-to-shallow transitions (GRB 050315, 050724), large X-ray flares (GRB 050502B, 
050724), fastly decUning (GRB 051210) and gradually declining (GRB 051221a, GRB 
050826; flux scale divided by 100 for clarity) afterglows. 



(Barthelmy et al. 2005b O'Brien et al. 2006). Three representative Swift X-ray 



lightcurves are shown in Figure |6] for both long and short GRBs. These X-ray 
hght curves start as early as a hundred seconds after the GRB trigger, and cover 
up to five decades in time. The complex behavior revealed in them sigificantly 
challenges traditional afterglow theoy, and calls into question some of the basic 
underlying assumptions. 

One of most striking results is that many of the early X-ray afterglows show 
a canonical behavior, where the light curve broadly consists of three distinct 
power law segments (Nousek et al. 2006). A bright rapid- falling where 
a > 3) afterglow immediately after the prompt emission dTaghaferri et al.|2005b| ) 
is followed by a steep-to-shallow transition, which is usually accompanied by a 
change in the spectrum power-law index. This is consistent with an interpretation 



(Nousek et al. 2006, Zhang et al. 2006) in which the first break occurs when the 



slowly decaying emission from the forward shock becomes dominant over the 
steeply decaying tail emission of the prompt 7-rays as seen from large angles 
( Kumar and PanaitescupOOO ) . Since these two components arise from physically 
distinct regions, their spectrum would generally be different. The shallow phase 
then transitions to the classical afterglow phase with no clear evidence for a 
spectral change. In some cases a jet break is seen at late times. The intermediate 
shallow flux stage is commonly interpreted as being caused by the continuous 



energy injection into the external shock (Nousek et al. 2006 Zhang et al. 2006) 



although orientation and complex jet structures have been also discussed as viable 
alternatives. The energy in the afterglow at these late times is estimated to be 
comparable to or smaller than that in the prompt gamma-ray emission, even when 
correcting for radiative losses from the afterglow shock at early times, implying a 
high efficiency of the prompt emission. The presence of the shallow decay phase 
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implies that most of the energy in the afterglow shock was either injected at late 
times after the prompt gamma-ray emission was over or was originally in slow 
material that would not have contributed to the prompt gamma-ray emission. 
This requires the prompt gamma-ray emission mechanism to be significantly more 
efficient than previous estimates. If a significant fraction of the radiated energy 
goes to photon energies above the observed range, the efficiency requirements of 
the prompt emission become even more severe. 



Table 1 Typical parameters of the canonical Swift X-ray light curve 









Approximate 


Phase 


Start T (s) 


Decay index" 


frequency 


Steep decline 


10^ - 10^ 


>3 


50% 


Shallow slope 


10^ - 10^ 


0.5 


60% 


Classical afterglow 


10^ - 10*^ 


1.3 


80% 


Jet break late phase 


10^ - 10*^ 


2.3 


20%^^ 


X-ray flares 


102 - 10^ 




50% 



"Decay index a defined by F = Fot 

^Of the 80% with no observed jet break, about half had afterglow obser- 
vations terminate before expected time of jet break. 



The average times, slopes, and frequencies characterizing these three distinct 
X-ray afterglow components are listed in Table [T} Most Swift X-ray light curves 
are broadly consistent with this basic temporal description, although in most 
cases we do not see all three power law segments, either because not all are present 
or because of limited temporal coverage. The large variety of behaviors exhibited 
by afterglows at different times in their evolution can be seen in Figure [7] which 
shows the temporal history for each individual afterglow as well as the evolution 
of the cumulative X-ray afterglow luminosity for a large sample of Swift events 
with known redshift. While broadly compatible with relativistic fireball models 
(Nousek et al. 2006 Zhang et al. 20061 the complex afterglow behavior that 
has been revealed poses new challenges of interpretation. The reader is refered to 



Granot (2008) for a more detailed account of the major strengths and weaknesses 



of the standard afterglow model, as well as some of the challenges that it faces 
in explaining recent data. 

Swift has also discovered flaring behavior appearing well after the prompt phase 
in ~50% of X-ray afterglows (Chincarini et al. 2007 Falcone et al. 2007 1 . An 
illustration of bursts with bright flares is shown in Figure [6} In some extreme 
cases, the late time flares have integrated energy similar to or exceeding the initial 
burst of 7-rays (Burrows et al. 2005b I . The rapid rise and decay, multiple flares 
in the same burst, and cases of fluence comparable to the prompt emission sug- 
gest that these flares are due to the same mechanism responsible for the prompt 
emission, which is usually attributed to the activity of the central powerhouse. 
When X-ray flares are observed by XRT, it is typically the case that no flaring is 



seen in the optical band by the UVOT. A notable example is GRB 060418 (Moli- 



nari et al. 2007 1, whose optical-infrared afterglow spectra is not consistent with 
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Figure 7 X-ray and optical lightcurves of GRB afterglows in the Swift era. (a) X-ray 
lightcurves of Swift burst afterglows. Data for long-duration bursts with known red- 
shifts, from GRB 050126 to GRB 070724A, have been gathered from the Swift XRT 



hghtcurve and spectral data depository at the UK Swift Science Data Centre (Evans 



et al. 2007 Evans et al.||2008 ). They are corrected uniformly to unabsorbed luminosity 



over 1.0 — 30.0 keV in the burst rest-frame, using the time-average afterglow spectrum, 
and plotted as a function of rest-frame time (pink lines) . Separately, afterglow lightcurve 
fits (Racusin et al. 2008b I, which exclude flaring intervals, are used to construct min- 
imum and maximum envelopes (black lines) and confidence intervals (gray bands) on 
the X-ray luminosities of the bursts as a function of rest-frame time: light gray regions 
delimit bands of 10% to 90% confidence, dark gray regions delimit bands of 25% to 75% 
confidence, and the median burst luminosity at any given time is shown (middle black 
line). (&) Optical lightcurves of Swift burst afterglows. Data for long-duration bursts 



with known redshifts and at least "bronze" quality published optical data (Kann et al. 



2007), from GRB 050408 to GRB070612A, are corrected uniformly to rest-frame f7-band 



luminosity using the inferred i?-band (z = 1) light-curves from Kann et al. (20071, and 
plotted as a function of rest-frame time (pink lines). Interpolated and "best fit" extrap- 
olated lightcurves are used to generate minimum and maximum envelopes and median 
luminosity estimates (black lines) and confidence intervals (gray regions), as in a. 
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a simple power-law extrapolation to soft X-ray energies and clearly requires two 
distinct spectral components. Not surprisingly, the broadband energy spectra of 
GRBs are complex and have spectra of wildly disparate shapes even for members 
belonging to a single class. 

Prior to Swift there were several reports of emission and absorption line ob- 
servations in the X-ray spectra of GRB afterglows. These included BeppoSAX 
observations of GRB 970508 ( |Piro et al.|[l999l ) and GRB 000214 dAntonelU et"al 
2000D , ASCA observations of GRB 970828 ( |Yoshida et ar||1999D , Chandra obser- 
vations of GRB 991216 ( |Piro et al.|2"000l ) and XMM observations of GRB 011211 
(Reeves et al. 2002 1 . None of the detections were of high statistical significance, 



but, combined, were of some credibility. Swift has not found any significant line 
features in comprehensive observations of the X-ray afterglow of more than 200 



GRBs (Romano et al. 20081. Hurkett et al. (20081 also made a very detailed 



Swift X-ray line search, with no positive result. These strong negative findings 
call into question the significance of previous results. 

A final, but widely discussed, observational development is the discovery in 
~25% of short bursts detected by Swift /BAT of an extended emission (EE) com- 
ponent lasting for ~100 seconds (Norris and Gehrels 2008). This component 



was clearly detected in HETE-2 burst GRB 050709 (Villasenor et al. 20051 and 



Swift burst GRB 050724 (Barthelmy et al. 2005c I . Archival searches have also 



found BATSE bursts with EE (Norris and Bonnell 2006 Lazzati et al. 2001 



Connaughton 2002). The EE is typically softer than the main peak and has an 
intensity range from 10^'^ to 10~^ times that of the initial short pulse complex. 
It is possible that many of the 75% of bursts without currently detected EE have 
this component at flux levels below detectability, although there are bursts with 
upper limits on the intensity of < 10~^ times that of the short pulse complex. 



3.2.2 Optical observations. With increasing frequency during the Swift 
era, optical observations of GRBs have been commencing almost immediately 
after 
to 



and in one noteworthy case (GRB 080319B; Racusin et al. 2008a I, prior 
the GRB trigger itself. These early optical observations have substantially 
enriched our appreciation for the complexity of the physical processes active 
during the prompt emission phase of GRBs, and during the afterglow that follows. 
The large diversity of optical afterglow lightcurves can be seen in Figure [7j 



Prompt emission. Observations prior to the Swift era demonstrated already 



that bright optical flashes such as that seen from GRB 990123 (Akerlof et al 



1999), were relatively rare (Kehoe et al. 2001 1. The success of the Swift mission 



has brought a vast increase in the rate of bursts accessible to rapid optical follow- 
up, and a corresponding increase in the number of events detected from early 
times, t < 100 s. The Swift UVOT itself routinely detects optical afterglows 
following the intial prompt slew of the satellite, with a 40% detection rate for 
such bursts ( Roming et al.|2009 ) that is only slightly lower than the 60% detection 
rate combining all observatories. 

The fastest routine responses to Swift alerts are realized by robotic ground- 
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based telescopes. The first discovery yielded by these observatories in the Swift 



era was of the 7-ray correlated component of the prompt optical emission (Ves- 



trand et al. 


2005 


Blake et al. 


2005 



not observed in every burst, but the mere fact of its correlation is sufficient to 
establish a common origin with the prompt 7-ray emission (e.g., internal shocks). 
When observed, the ratio of the correlated 7-ray to optical flux densities has been 
found to be roug hly 10^ to one. 

In contrast to bursts with 7-ray correlated emission, the common burst is 
now revealed to either exhibit a single power-law decay from early times (Rykoff 

or to exhibit a flat or 



et al. 


2005b Quimby et al. 


2006a 




Yost et al. 


2006a 


1, 


rising 


(Rykoff et al. 


2004a 


] 


lykoff et al. 


2006a 1 


or rebri 



2007) optical light curve before it enters the standard power-law afterglow decay. 
The initial brightness of the typical counterpart is V ~ 14 to 17 maa^ (Roming 



et al.|[2009 |, which has made observations challenging for smaller (<lm) robotic 



facilities, and has limited the extent of the lightcurves collected by the UVOT. 
A few new observations of bright flaring optical emission have been collected. 



which are usually interpreted as emission from the reverse shock region (Sari 



and Piraii||1999 Meszaros and Rees||1999 ). The early optical/NIR emission from 



GRB 041219 would have rivalled that seen from GRB 990123 if not for the large 
Galactic extinction along the line of sight (Vestrand et al. 20051. Of the three 
distinct peaks observed by PAIRITEL, the second may represent the onset of 
the afterglow (reverse shock) contribution ( Blake et al.||2005 ) - if so, a relatively 
small Lorentz factor, T ~ 70, is derived by associating the flaring peak time with 
deceleration time of the relativistic blast wave. Observations of GRB 050525A 
with UVOT dBlustin et ar|[2006l ) and GRB 06011 IB with TAROT (in a unique 
time-resolved tracking mode; Klotz et al.| 2006 ) show the "flattening" light-curve 
famihar from GRB 021211 ( [Zhang et al]2003[ ). Intriguingly, the "high redshift" 



GRB 050904 ( [TagUaferri ^ral]|2005a| [Cusumano eTaL][2006l ) at z = 6.29 ( [Kawai 



et al.| 2006 1 also had prompt optical emission observed ( Boer et aL]|2006 Haislip 



et al. 20061, with a brightness, single-pulse structure, and fast-fading behavior 



reminiscent of GRB 990123 and thus, potentially also interpreted as reverse shock 
emission. 

Swift detection of the "naked eye burst" GRB 080319B, which peaked at vi- 
sual magnitude V = 5.3 has now delivered the richest dataset, by far, addressing 
the prompt optical emission and its evolution into a standard fading afterglow 
dRacusin et "aL]|2008a[ [Bloom et aL][2009l jWozniak et"al][2008| . This is only par- 
tially due to the extreme brightness of the event; the fact that it occurred at an 
equatorial location, in the night sky above the Western hemisphere, is probably 
even more important. The fact that it occurred within just one hour and ten 
degrees of the preceding GRB080319A means that it is also the only event with 



strong contraints on optical precursor emission from pointed telescopes (Racusin 



et al. 2008a|. The GRB080319B dataset is rich enough that its ramifications 
are still being grappled with. In an overall sense the prompt optical emission 
correlates well with the 7-ray lightcurve; in detail, though, the individual pulses 



The ~ 60% of those counterparts not detected are presumably fainter than this. 
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observed in both bands do not track precisely, probably suggesting the presence 
of at least two distinct cooling processes at play within the dissipation region 
( Racusin et al.|2008a |.GRB 0803 19B also allowed astronomers to rule out inverse 
comptonization as a relevant mechanism (Piran et al.||2009 Zou et ar]|2009 1 . 



Afterglow emission. The observed properties of the optical afterglow are 



largely familiar from observations prior to the Swift era (van Paradijs et al. 112000 



Meszaros 20021 and will not be reviewed in detail here. An intriguing feature 



of later-time afterglows, revealed by the rich multiband lightcurves available in 
the Swift era, has been the occasional presence of chromatic lightcurve breaks, 
where the X-rays show a clear break (steepening of the flux decay rate) while 
the optical does not. The break in the X-ray light curve is usually identified 
with the end of the shallow decay phase. The optical lightcurve follows a sin- 
gle power law decay, usually with a temporal decay index intermediate between 
those in the X-rays before and after the break. Accommodating such chromatic 
breaks in a model where the X-ray and optical emissions arise from the same 
emitting region (Panaitescu et al. 20061 requires not only a temporal evolution 
of the underlying microphysical conditions within the emitting region but also 
fine tuning their photon arrival times in such a way that a break in the X-ray 
will be produced but not at optical wavelengths. Alternatively, the X-ray and 
optical photons may arise from physically distinct regions, which would naturally 
account for their seemingly decoupled behavior. Observations of the naked-eye 
GRB 080319B have sharpened this debate. This is because the optical and X-ray 
decays in this event exhibit discrepant behavior over two orders of magnitude in 
time, from 100s to a few times 10^ s post-burst ( Racusin et al.|2008a Kumar and 



Panaitescu 20081. 



Recent results showed that during the initial 500s of observations 15% of UVOT 
lightcurves are seen to rise - with an average peak time of 400s, 58% decay from 
the onset of observations and the remainder are consistent with being flat, but 
could be rising or decaying (Oates et al. 2009). This leads to a wide range of 
temporal indices measured before 500s, —1.17 < a < 0.21. Such behavior is also 



observed by ground based telescopes (Rykoff et al. 2005a Quimby et al. 2006b 



Yost et al. 2006b, Rykoff et al. 2006b) and was seen in pre-Swift observations 



(Rykoff et al. 2004b). No color evolution was observed during the rising phase 
of the UVOT lightcurves. One likely scenario is that the rise is caused by the 
jet plowing into the external medium - the start of the forward shock. If this 
scenario is correct, then it is possible to determine the Lorentz factor of the jet at 
the time of the peak, giving T ~ 180 and a lower limit T > 230 for those without 
observed peaks (Oates et al. 2009). 

After a few hundred seconds the afterglow decays as a single power-law with a 
temporal index, measured after 500s, of —1.20 < a < —0.52. This is consistent 
with pre-Swift observations and the range is similar to the shallow decay of the 
XRT canonical model. However, 20% of UVOT lightcurves are seen with a broken 
power-law decay after 500s. For these afterglows, the first decay has a range 



21 



between —0.74 < a < —0.46, which is most consistent with the shahow decay 
of the X-ray canonical hghtcurve and the second decay has a temporal range of 
—1.72 < a < —1.34, which is most consitent with the classic afterglow phase. 

At late times, t ^ 1 day, it was common in the pve-Swift era to observe a 
steepening of the optical decay to power-law indices ao > 2, and to associate this 
epoch with the "jet break" transition of the underlying, relativistically expanding 
jet. As jet breaks are intrinsically achromatic, they are expected to manifest in 
the X-ray light curves from Swift as well as in ground-based optical observations. 
In fact, over the first three years of Swift operations only a handful of convincing 



jet-break candidates were identified ( 


Blustin et al. 


2006 


Stanek et al. 


2007 


Dai 


et al.|2007 


Willingale et al.||2007 


Kocevski and Butler|2008 


) , leading to concerns 



about the viability of this picture. Some bursts with very long coverage have been 



convincingly shown to possess no significant breaks, e.g., GRB 060729 (Grupe 



etal][2007| ). 

Since then, however, deep optical imaging observations have revealed evidence 
for jet breaks in several additional Swift afterglows, a significant fraction of those 



monitored to late times (Dai et al. 20081. Systematic analysis of the Swift XRT 



data has revealed strong evidence for a jet-break signature in 12% of the X-ray 
lightcurves ( Racusin et al.||2008b ). A consensus has developed that jet breaks for 
typical Swift bursts may be occurring at late times and faint flux levels that are 
beyond the limit of the standard ground- and space-based campaigns. Separately, 
it has been suggested that a distinct additional spectral component could hide 
the jet-break signature in the X-ray band. For example, in the case of the bright 



GRB 070125 (Chandra et al. 2008), inverse Compton scattering of the synchrotron 
optical photons has been put forward as an explanation for the missing jet-break 
feature. 



3.2.3 Radio observations. Radio afterglow observations are unique in hav- 
ing led to both indirect and direct demonstrations of relativistic expansion, via 



scintillation (e.g. Chandra et al. 20081 and VLBI (Taylor et al. 20051 observa- 
tions. They provide access to the properties of bright afterglows on the smallest 
angular scales. The limited number of sensitive high-resolution radio facilities, 
and the sensitivity limits of those facilities, have prevented a proportional ex- 
ploitation of the greatly-increased burst rate from Swift, primarily because the 
compensating feature of this increased burst rate has been a greater median red- 
shift and a lower characteristic afterglow flux. 

Radio afterglow observations nonetheless continue to play a vital role in accu- 
rately estimating blastwave kinetic energies dOren et al.|2004||Granot et al.|2005a| 



Kaneko et al. 20071, with radio detections contributing crucially to the demon- 

-j^q52 gj,g^ kinetic energy associated with the 



stration of the extremely large {E 
high-redshift GRB 050904 ( [Frail et al.|[2006l |Gou et al.|[2007l ), and to constrain- 
ing the relativistic energy associated with the nearby GRB060218/SN2006aj 
( Soderberg et al.||2006b[ ). Radio data can also provide a crucial "third check" on 
claims of jet-break detections, as with the broadband afterglow models applied 
to GRB 050820 ( |Cenko et al.][2006| ) and GRB 070125 ( [Chandra et al][2008l ). 
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Figure 8 The growth of GRB 030329 with time as measured using VLBI (by Pihlstrom 
at al.|2007 and references therein) . In the background are the images from (a) April 2003 
(15 GHz), (b) June 2003 (8 GHz), (c) November 2003 (8 GHz), and (d) June 2005 (5 
GHz) with the same intensity scale. The resolution for the four images is not constant 
in time, but is accounted for in the analysis of the source size. 



For relatively nearby GRBs which may be associated with a supernova, radio 
observations have proven invaluable in providing evidence of the physical expan- 
sion of GRB/SNe ejecta through direct imaging. Figure [s] shows the expansion 
from ~ 0.02 mas to ~ 0.35 mas of the radio image of GRB 030329 over a time 
span from ~ 15s after the GRB trigger to ~ lO^s (Pihlstrom et al. 2007 1 . 



3.3 Interpreting Prompt and Afterglow Emission 

The isotropic-equivalent luminosity of GRB X-ray afterglows scaled to t 11 
hours after the burst in the source frame, can be used as an approximate estima- 



tor for the energy in the afterglow shock for the following reasons (Freedman and 



Waxman 


2001 


Piran et al. 


2001 



band is typically above the two characteristic synchrotron frequencies, so that the 
flux has very weak dependence on microphysical parameters and no dependence 
on the external density, both of which are associated with relatively large uncer- 
tainties. Second, at 11 hr the Lorentz factor of the afterglow shock is sufficiently 
small (r ~ 10) so that a large fraction of the jet is visible (out to an angle of 
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Figure 9 Isotropic-equivalent luminosity of GRB X-ray afterglows scaled to t — 11 hr 
(5 keV source frame) after the burst trigger as a function of their isotropic 7-ray energy 
release (adapted from Nyse wander et al.]|2008 1 . 



~ ~ 0.1 rad around the line of sight) and local inhomogeneities on small 
angular scales are averaged out. Finally, the fact that the ratio of Lx(ll hr) and 
-Eiso is fairly constant for most GRBs, suggests that both can serve as a reasonable 
measure of the isotropic-equivalent energy content of the ejected outflow. 

Figure [9] shows Lx(ll hr) at 5 keV rest-frame energy as a function of their 
isotropic 7-ray energy release for a large sample of GRBs. A linear relation, 
Lx oc (11 hr) oc E^^iso, seems to be broadly consistent with the data, probably 
suggesting a roughly universal efficiency for converting kinetic energy into 7— rays 



in the prompt emission for both long and short GRBs ( 


Lee et al. 


2005a 


Nousek 


et al.||2006 


Bloom et al.||2007 


Kaneko et al.||2007 


Nysewander et al.||2008|. This 


"universal" efficiency is also likely to be high (i.e., the remaining kinetic energy 



is comparable to, or even smaller than, the energy dissipated and radiated in the 
prompt emission). If this is the case, the well-known efficiency problem for long 
GRBs also persists for short events. 



3.4 The Radiated Energy Inventory 

With the advent of Swift, the observational inventory for GRBs has become rich 
enough to allow estimates of their radiated energy content. A compilation of the 



radiated energy in both the prompt and afterglow phases is presented in Figure 10 
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Figure 10 The radiated-energy inventory of Swift GRBs. Left Panel: Summary of the 
isotropic-equivalent total emitted energy o f the pro mpt and afterglow emission for Swift 
GRBs in the source frame (adapted from Kaneko et al. 2007). Right Panel Bottom: Dis- 
tributions of cumulative isotropic-equivalent total emitted energy for Swift GRBs. Right 
Panel Top: Comparisons of coUimation-corrected total emitted 7-ray energy, i?iso(fi/47r), 
of pre-Swift GRBs, where Eiso is used as an upper limit for GRBs with no jet angle con- 
straints. 



To investigate the energy dissipation behavior in the X-ray and optical afterglow, 
we fitted a natural cubic spline function to the afterglow histories (shown in 
Figure [7]) for each individual afterglow and estimated the cumulative emitted 
energy. The start and end times of the integration were the first and the last 
points of the actual observations. The result is a global portrait of the effects of 
the physical processes responsible for GRB evolution, operating on scales ranging 
from AU to parsec lengths. The compilation also offers a way to assess how well 
we understand the physics of GRBs, by the degree of consistency among related 
entries. 



Figure 10 shows that the isotropic X-ray emission, Ex.iso, for most Swift GRBs 
spans the range of ~ 10^*^ to 10^^ erg, which is comparable to that emitted 
during the prompt 7-ray phase. Not surprisingly, events that have large isotropic- 
equivalent energy in 7-rays have larger i?x,iso) indicating a reasonably narrow 
spread in the efficiency of converting the afterglow kinetic energy into radiation. 
As can be seen also in Figure 10, the isotropic equivalent energy that is radiated 
at optical wavelengths is ~ 2 orders of magnitude smaller than that in X-rays 
and 7-rays. This is predominantly due to the fact that vFi^ typically peaks closer 
to the X-rays than to the optical, and it is very flat above its peak while it falls 
much faster toward lower energies. Finally, since these are isotropic equivalent 
energies, most of the contribution to the radiated afterglow energy, especially 
at radio wavelengths, is from significantly later times than for Exiso, and the 
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collimation of the outflow together with relativistic beaming effects could result 
in much larger -Ex.iso than -Er^so- 



Of the three nearby SN-GRB events plotted in Figure [10) only GRB 030329 
falls within the ~ lO^'^ erg range, while the other two events fall between ~ 10^^ 
and ~ 10^^ erg (Kaneko et al. 2007 1 . We note that there exists a selection effect 



based on the observed photon flux: an event is more likely to be detected when 
it is closer to us than farther, for a given intrinsic luminosity. 

One of the liveliest debated issues associated with GRBs concerns the total 
energy released during the explosion: are GRBs standard candles? The GRB 
community has vacillated between initial claims that the GRB intrinsic luminos- 
ity distribution was very narrow ( Horack et al -111994 1 to discounting all standard- 



candle claims, to accepting a standard total GRB energy of ~ 10 ergs (Frail 



et al. 



20011, and to diversifying GRBs into normal and sub-energetic classes 



(Soderberg et al. 2008, Ramirez-Ruiz et al. 2005a I . The important new devel 



opment is that we now have significant observational support for the existence 
of a sub-energetic population based on the different amounts of relativistic en- 
ergy released during the initial explosion. A network of theoretical tests lends 



credence to this idea (Waxman 2004a 


Kaneko et al.||2007 


Granot and Ramirez- 


Ruiz 


2004 


Granot and Kumar||2006 




Plan et al. 


2006 


). The existence of a wide 



range of intrinsic energies could further challenge the use of GRBs as standard 
candles. 



4 ENVIRONMENTS AND HOST GALAXIES 

Much of what we know about GRBs has been derived not from observations of 
the prompt burst radiation itself, but from studies of their afterglows - as they 
illuminate the circumburst surroundings - and their host galaxies. In the sections 
that follow we discuss the primary insights that have been derived in this manner 
from the study of GRB environments. 

4.1 Cosmological Setting 

In Figure [4] we present the redshift distribution of all Swift-detected gamma-ray 
bursts. Swift and other current missions observe GRBs to cosmological distances 
quite readily; indeed, the three highest-fluence, known-redshift bursts observed 
by Swift have been at z = 0.61 (GRB 050525A), z = 2.82 (GRB 050603), and 
z = 1.26 (GRB 061007) - already spanning 40% of cosmic history. Historically, 
the majority of redshifts have been collected via host galaxy spectroscopy; in the 
Swift era, however, this pattern has been reversed - except for the short bursts 
- with the great majority of redshifts now being derived from afterglow spectra. 



In addition to theoretical arguments that posit different physical origins for 



short and long burs 
associated supernov; 


ts (e.g., Katz and Canel 


1996|), the absence of short burst- 


le to deep limits (|Hjorth et al.||2005a 


Fox et al.||2005 


Bloom 


and Prochaska 2006 


Castro- Tirado et al. 2005 


1, 100 times fainter than SN 1998bw 



in the best cases, argues for a distinct origin of the short and long bursts. In 
agreement with this picture, the redshift distributions of the two populations are 
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not consistent. 

For the long bursts, which are associated with active star formation and, in 
particular, the deaths of massive stars, it is interesting to explore whether their 
distribution in redshift is consistent with other measures of cosmic star forma- 
tion. The greatly-increased number of redshifts available for Swift bursts have 
motivated several such comparisons, and for the first time, estimates of cosmic 



star formation at high redshift z > 4 using the Swift redshift sample ( Chary et al. 



2007, Yiiksel et al. 20081. 



The star-formation studies over 1 < z < 4 confirm, in a broad sense, that 
the GRB redshift distribution remains consistent with independent measures of 
star formation (Jakobsson et al. 2006). However, there are signs of differential 



evolution of the GRB rate, in the sense that the GRB rate increases more rapidly 
with increasing redshift than expected based on star-formation measures alone 



(Le and Dermer 2007 Guetta and Piran 2007 Kistler et al. 2008). This evi 



dence, currently present at roughly 95%-confidence, may strengthen significantly 
in coming years. If so, this would provide a sign of bias towards low-mass and 
low-met allicity host galaxies - and potentially, low-metallicity progenitors - for 
the long-duration bursts (Section |4.2[ ). The theoretical curves accompanying the 
GRB redshift distribution (Figure |4]) which show the evolution of a comoving 
volume element in the universe, and the volume convolved with star formation 
rate, appear to indicate that the observed distribution is wider than expected. 

The short-burst redshift distribution, so far drawn exclusively from host-galaxy 
observations, has also been compared to star-formation metrics. In this case, 
however, the intent has been to explore "time-delayed" progenitor models that 
correlate with star formation through a parametrized (log-normal or power-law) 
delay function ( Nakar et al.|2006a Guetta and Piran|2006 Salvaterra et al.|2608 ). 
Consistent with the relatively large fraction of events at z ^ 0.5 compared to long 
bursts, and with host galaxy demographics ( Zheng and Ramirez-Ruiz||2007 Shin| 
and Berger 20071, these studies have concluded that a long-lived (r ^ 1 Gyr) 



progenitor is required for these models to be consistent with Swift-era redshift 
measurements and the distribution of short burst fluences from BATSE. This in 
turn has led to relatively high estimates of the volumetric local short burst rate. 



at least an order of magnitude greater than the local rate of long bursts (Guetta 



et al. 2005), and a correspondingly optimistic set of predictions for Enhanced 



LIGO, VIRGO, and other ground-based gravity- wave detectors. 



4.2 Host Galaxies of Long Bursts 



Surveys of GRB host galaxies in the pre- Swift era ( Le Floc'h et al.|[2003 Fynbo 
et al. 2003 ) necessarily focused on the host galaxies of long-duration bursts. These 
surveys established a standard picture for the GRB hosts as sub-L* galaxies 
(median L ~ O.IL*) with exponential-disk light profiles (Conselice et al. 2005 



Wainwright et al.]|2007 ) and high specific star-formation rates (SSFR ~ 1 Gyr ^) 
(Christensen et al. 2004). A selection of GRB host galaxies, as imaged by HST, 



are shown in Figure 11 



This picture has not substantially changed in the Swift era. To the contrary. 
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Figure 1 1 A selection of the host galaxies of long-duration (top row) and short-duration 
(bottom row) gamma-ray bursts, as imaged by the Hubble Space Telescope. An attempt 
has been made to choose pairs of long and short burst host galaxies with comparable 
redshifts; lower-redshift hosts are emphasized as these reveal their structure more readily 
in short exposures. Images are oriented with North up and East to the left, and the 
physical length scale for a one-arcsecond angular distance is indicated in each panel 
(except for GRB 051227); arrows point to the location of the burst where this is known 
to ~pixel precision. Individual burst notes: GRB 030329 was the first classical long GRB 
to be associated with a well-observed spectroscopic supernova ( Stanek et al.|200"3 Hjorth 
et al.|2003| ); XRF 020903 was the first X-ray Flash event to yield a redshift measurement 
Soderberg et al. ( 2004 ) ; GRB 050709 was the first short burst with optical afterglow - 



indicated by the arrow - detected Fox et al. (20051, Hjorth et al. (2005b I; GRB 050509B 



was the first short burst with detected afterglow ( Gehrels et al.|2005 Bloom et al.||2006 l 
GRB 051227 has a faint candidate host, of unknown redshift probably greater than 1, 
visible at the optical afterglow location; the spiral galaxy to the east has redshift z = 0.714 



Foley et al. (2005). Long-burst host images from Wainwright et al. (20071; short- burst 



host images from Fox et al. ( 2005 1 and this work. 



the result of various ground- and space-based efforts to characterize GRB host 



galaxies ( 


Le Floc'h et al. 


2006 


Chary et al. 


2007 


Fynbo et al. 


2008 


Savaglio 


et al. 


20091 have confirmed this basic outline and expanded its domain of ap- 



plicability to high redshift, combining the power of the Swift burst catalog with 
Spitzer observations. At higher redshifts, z ^ 3, it seems particularly interest- 
ing to use the GRB host galaxies to explore the evolution of mass-met allicity 
relationships that are typically compiled using field galaxies at low redshift and 



high- mass galaxy samples at high-redshift (Figure 12). Since GRB host red- 
shifts are typically secured via afterglow spectroscopy, the hosts themselves are 
uniquely free of mass and luminosity selection effects. In this context, even up- 
per limits provide useful constraints on mass-met allicity correlations and their 



evolution with redshift (Chary et al. 2007 Savagho et al. 2009 Berger 2009). 



Such studies have also served to place GRB host galaxies in the context of other 
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Figure 12 Metallicity as a function of i3-band absolute magnitude for the host galaxies 
of short (red) and long (blue) GRBs. The yellow bars mark the 14 — 86 percentil e range 
for galaxies at z ^ 0.1 from the Sloan Digital Sk y Survey (Tremonti et al.|2004 l, while 
crosses designate field galaxies at z ~ 0.3 — 1 ([Kobulnicky and Kewley |2004p . Both 
field samples exhibit a clear luminosity-met allicity r elatio n. I 'he long GKB hosts tend 
to exhibit lower than expected metallicities ( Stanek et aH2006[ ), while the hosts of short 
GRBs have higher metallicities by about 0.6 de x and are nior eover in excellent agreement 
with the luminosity-metallicity relation. From Berger (20091. 



high-redshift galaxy populations incuding the Lyman-break and Lyman-alpha 
emitting galaxies ( Fynbo et al.||2008 |. 

Finally, it is somewhat reassuring that numerical simulations of star-forming 
galaxy populations can generate galaxy subpopulations that reproduce the basic 



properties of the GRB host galaxies in several important respects (Courty et al. 



2004, Courty et al.||2007 ). In particular, selecting for high specific star- formation 



rate generates mock galaxy catalogs with similar masses, luminosities, and colors 
to GRB host galaxies. 



4.2.1 Metallicity matters. While GRB host galaxies are often studied for 
the insights they provide about larger astrophysical questions, including the his- 
tory of star formation through cosmic time, they also have the potential to shed 
light on the nature of the GRB progenitors. The association of long-duration 
bursts with star formation, for example, was proposed after observation of just 
two host galaxies ( Paczyhski|1998 1 , and demonstrated firmly from the properties 
of the first twenty (Bloom et al. 2002). 



Recent years have seen a surge of interest in the question of whether GRB host 
galaxies, and hence, presumably, GRB progenitors, are metal-poor by comparison 
to the larger population of star-forming galaxies. This question has been explored 
from a variety of perspectives, and the result of these studies, still under active 
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debate, may eventually help refine our picture of the massive stellar death behind 
each GRB. 

The metallicities of GRB host galaxies, and indeed, detailed abundance profiles, 
can be measured directly from high-resolution spectroscopy of bright afterglows. 
This area has seen dramatic progress in the Swift era, with prompt arcsecond 
positions and visual magnitude measurements from Swift UVOT and ground- 
based robotic telescopes feeding rapid-response spectroscopy from large-aperture 
facilities. One chief result of these efforts have been the collection of detailed 
abundance characterizations for multiple bursts ( Vreeswijk et al.|2004 Chen et al. 



20051 IFynbo et al ."2006at [Bloom et ar|[2007l [Prochaska eral]|2007l [Chen et al. 



2007| [Price et al.||2007, ,Chen et al.[|2008[ ). In addition, the metallicities of GRB 



host galaxies have been measured via standard emission-line diagnostics ( Stanek 



et al. 2006 Thone et al. 2008 SavagHo et al. 20091, where possible, and in some 



cases tentative conclusions have been drawn on the basis of cruder relations such 



as mass-met allicity metrics (Berger et al. 2007a I and host galaxy luminosities 
and morphologies ( Fruchter et al.|[2006 |. 

The conclusions of these studies have yet to be reconciled into a single coherent 
picture of the nature of the GRB host galaxies and their relationship to other 
low- and high-redshift galaxy populations. However, the wealth of data do serve 
to define some associated issues with reasonable clarity. First, the metallicities of 
GRB host galaxies at z ^ 1 are significantly {Z ~ O.IZq) sub-solar (Savaglio et al 



20091, consistent with the sub-solar metallicities measured for GRB host galaxies 



via absorption spectroscopy at z ^ 2 (e.g., Chen et al. 2008). These sub-solar 
metallicities are neither surprising nor unusual for galaxy populations at high 
redshift (Fynbo et al. 2008 Savaglio et al. 20091; moreover, several candidate 



higher- metallicity hosts have been identified, although not yet confirmed (Berger 



et al.[[2007a[ [Fruchter et 'al][2006[ ) . At the same time, GRB host galaxies seem to 



be readily distinguished, in luminosity and morphology, from the host galaxies of 
core-collapse supernovae at similar redshifts (Fruchter et al. 20061, and the host 
galaxies of the lowest-redshift z ^ 0.2 bursts have uniformly low metallicities that 
strongly distinguish them from the bulk of the low-redshift galaxy population 
(Stanek et al. 2006 Thone et al. 2008), and indeed, from the host galaxies of 
nearby type Ibc supernovae (Modjaz et al. 2008). 

These somewhat divergent findings might be reconciled in a picture where the 
GRB progenitors prefer (or require) a low-met allicity environment, since the in- 
creasing prevalence of such environments at z ^ 1 would allow GRB host galaxies 
to present an increasingly fair sample of the population of star-forming galaxies 
at these higher redshifts. This argument would also dovetail with observations 
that the GRB rate seems to increase with redshift faster than the cosmic star 



formation rate, as mentioned above (Le and Dermer 2007 Guetta and Piran 



2007 Kistler et al. 2008). However, the claim that GRB host galaxies represent 



a fair sample of star-forming galaxies, even at z 1, remains in dispute (e.g.. 
Figure [l2j |Berger[[2009| ) . 

A possibly significant implication of a metallicity-dependent GRB rate would 
be an offset between the true star-formation rate and that traced by GRBs. 
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If GRBs in low-metallicity environments and low-mass galaxies are more lumi- 
nous, then they are likely over-represented in GRB samples. Low-mass galaxies 
and galaxy outskirts have lower metallicity on average and thus may yield more 



(and/or more luminous) GRBs compared to high-mass galaxies (Ramirez-Ruiz 
et al. 2002b| ). As galaxy mass builds up through mergers, it is also possible that 
the highest-z GRBs could be systematically more luminous due to their lower- 
mass host galaxies, an intriguing hypothesis given the extreme luminosities of 
some of the highest-redshift bursts of the Swift era, including GRB 050904. 

Finally, it is worth noting that high-resolution afterglow spectra show absorp- 
tion features imprinted on the afterglow by gas at multiple, widely divergent phys- 
ical scales, possibly extending from d ~ 10 pc to Gpc distances (Section [4.2.2 ). 
With the GRB metallicity question ultimately referring to the nature of the pro- 
genitor itself, searches for definitive signatures of the progenitor's stellar wind 



material ( 


Moller et al.| 


2002 


Schaefer et al. 


|2003 


Mirabal et al. 


2003 


Bloom 


et al.||2007 


Chen et al. 


2007 


Fox et al. 


2008 


Prochaska et al. 


2008 


1 should con- 



tinue to be pursued and refined. At the same time, in discussing the host galaxies 
as a population, more common galaxy-integrated, emission-line diagnostics may 
better serve to place GRB hosts in the proper cosmological context. 



4.2.2 Subgalactic environments. The same fast-response spectra that 
have enabled characterization of elemental abundances in GRB host galaxies have 
led to a series of discoveries regarding absorbing gas structures on sub-galactic 
scales within the GRB host: a rich array of high-ionization fine-structure Fe tran- 
sitions in GRB 051111 ( Berger et al.|2005a Penprase et al.|2006 Prochaska et al 
2006); discovery of time- variability of such fine structure features in three bursts 



demonstrating excitation by UV photons from the burst flash and young after- 
glow ( [Delia Valle et al.|2006b[ [Vreeswijk et al.|2007[ [D'Elia et al.|2008[ ) ; discovery 
of high-ionization N V features, providing evidence for absorption by gas within 
d <^ lOpc of the burst ( Fox et al.|2008 Prochaska et al.[[2008 ); and most recently, 
the first detection of molecular gas along the line of sight to a GRB afterglow 



(Figure [Tsj [Prochaska et al.[[2009[ ) . 

The range of these discoveries reveals a surprising complication in the inter- 
pretation of afterglow spectra. Depending on the burst, it may be necessary to 
account for absorbing structures on every scale from the immediate parsec-size 
circumburst environment, to surrounding or intervening molecular clouds, to host 
galaxy absorbers on kiloparsec scales, and ultimately, to low-ionization gas associ- 
ated with cosmological structures at gigaparsec distances. Similar complications 
are a feature of quasar absorption line studies; there, however, the situation is 
simplified by the persistent nature and overwhelming ionizing power of the quasar 
itself. 

Indeed, the situation may be even more complicated than this. The discovery 
that the frequency {dN/ dz) of strong Mg II absorbers along GRB lines of sight 



seems to be roughly four times the frequency along quasar lines of sight ( Prochter 



et al. 2006) suggests that some or most such absorbers seen toward GRBs - 
despite their large separation in redshift from the host galaxy - may be intrinsic 
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Figure 13 Keck/LRIS spectrum of the afterglow for GRB 080607 (Prochaska at al. 20091 



The red dashed hnes indicate a model of the intrinsic afterglow spectrum reddened heavily 
by dust in the host galaxy (rest-frame Ay « 3.2). At A « 4900 A one identifies a damped 
Lya profile associated with H I gas near the GRB. The shaded region overplotted on the 
data corresponds to an H I column density = 10^^'^*°'^^cm— 2. The model (green 
solid line) includes absorption from H2 Lyman- Werner transitions. The line opacity at 
A > 5500 A is dominated by metal-line transitions from gas in the host galaxy and 
includes bandheads of the CO molecule. Suprisingly, this is the only sightline to date 
to show strong molecular absorption ( Tumlinson et al.]|2007 |. It also exhibits a roughly 
solar metallicity. The figure also shows the spectral region or features corresponding to 
intergalactic Lya and Mgll absorption. The redshift for GRB 080607 is 2; = 3.036. 



to the GRB environment ( Porciani et aL|[2007 |. 



4.3 Host Galaxies of Short Bursts 



At the time of the Swift launch, the greatest mystery in GRB astronomy was 
the origin of short GRBs. A major step forward was made in summer 2005 with 
the localization and afterglow detection of three short bursts, GRB 050509B, 
050709 and 070724. These events were found to be localized in regions of low 
star formation, either in low star forming elliptical galaxies as for GRB 050509B 
dGehrels et al.|2005| [Bloom et al.|2006l ) and GRB 050724 ( [Barthelmy et al.|2005c 



et al. 



2005 



Berger et al.|2005b | or in a region of a galaxy with low star formation (IVillasenor 



Fox et al. 2005 Hjorth et al. 2005b I . This was in stark contrast to 



long bursts which are associated with star forming regions, and implied a non 



32 



massive core-collapse origin. 

More than three years after these first short burst localizations, the catalog 
of confidently-identified short burst host galaxies is growing to the point where 
systematic studies can be carried out (Berger 2009) - although the effects of 
uncertain burst attribution (i.e., long or short?), uncertain host identification 
(especially for bursts with only Swift BAT or XRT localizations), and unknown 
redshifts for faint candidate hosts conspire to keep any conclusions largely quali- 
tative at this time. 

Indeed, without direct afterglow spectroscopy, association of short bursts with 
candidate host galaxies and host galaxy clusters must be approached probabilis- 
tically. In cases where a well-localized (preferably sub-arcsec) afterglow falls on a 
luminous region of the candidate host, or within a high- mass or high-redshift clus- 
ter, the association can probably be considered secure; however, in other cases an 



a posteriori estimate of the probability of association must be made ( Bloom and 



Prochaska 2006 Fox and Roming 2007). Such estimates are inevitably strongly 



dependent on input assumptions. For example: What lifetime and kick velocity 
distributions might be appropriate for progenitor binary systems? What about 
other possible progenitor classes? Any assumptions must be carefully considered 
before and after they are applied. 

With these caveats, a picture of the short burst host galaxy population as a 
whole has developed (Figure 11 ). It consists of three classes of host, two of which 
became apparent soon after the short burst afterglow revolution of 2005: the 
low-redshift (z ^ 0.5), high-mass (L ~ L*), early-type host galaxies and galaxy 

|2005[|Fox et al.||2005[ [Bloom et al.||2006 



clusters, on the one hand (Gehrels et al. 



Berger et al.]|2005b| |Gal-Yam et "aL][2008D , typified by the hosts of GRB 050509B, 
GRB 050724, and GRB 050813; and the low-redshift, sub-L*, late-type galaxies 
on the other, typified by the hosts of GRB 050709 ( |Fox et al.|[2005|Covino et al. 



2006| ), GRB 051221 [Soderberg etal.l ( |2006a| , and GRB 061006 [D'Avanzo etaL 



( |2009D . 

The third class of short burst host galaxies consists of faint, star- forming galax- 



ies at z ;^ 1 (Berger et al. 



2007b 



Cenko et al. 



2008 



Berger 2009), reminiscent 



of the host galaxies of long bursts. The existence of such higher-redshift, star- 
forming short burst host galaxies was predicted ( Belczynski et alT]|2006 ) based on 
binary population synthesis models assuming the compact object merger model 
for the short bursts. These simulations yield a bimodal distribution of lifetimes 
for merging systems, with a spike of mergers at short timescales, r ^ 100 Myr, 
followed by a dominant merger population with a lifetime distribution. 

The association of short GRBs with both star-forming and early-type galax- 
ies has led to analogies with type la supernovae, whose host demographics have 
similarly provided evidence for a wide distribution of delay times between for- 
mation and explosion. At the same time, both core-collapse supernovae and 
long-duration GRBs are observed (almost) exclusively in late-type star-forming 
galaxies. As with supernovae and long-duration bursts, a detailed census of short 
burst redshifts, host galaxy types, and burst locations within those hosts will 
undoubtedly help to constrain progenitor models. 
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4.4 Neither Long Nor Short 



Some interesting gamma-ray burst host galaxies are not obviously associated with 
either the long or short burst classes. In fact, the "peculiar" cases of GRB 060505 
and GRB 060614 probably provide the first examples where studies of the host 
galaxy properties have been applied to argue for a (long- or short-) nature of the 



bursts themselves ( [Fynbo et al.||2006b| |Gal-Yam et al.||2006[ |Ofek et al.||2007 
Thone et al.]|2008[ ). 



GRB 060614 was a low-redshift, long-duration burst with no detection of a 
coincident supernova to deep limits. It was a bright burst (fluence in 15 — 150 
keV band of 2.2 x 10~^ erg cm~^) and well-studied in the X-ray and optical. 
With a T90 duration of 102 s, it seems to fall squarely in the long burst category. 



A host galaxy was found (Gal- Yam et al. 2006 Fynbo et al. 2006b, Delia Valle 



et al. 2006a I at z = 0.125 and deep searches were made for a coincident supernova. 



All other well-observed nearby GRBs have had supernovae detected, but this one 
did not to limits > 100 times fainter than previous detections. 

GRB 060614 shares some characteristics with short bursts ( Gehrels et al.|2006 |. 
The BAT light curve shows a first short, hard-spectrum episode of emission (last- 
ing 5 s) followed by an extended and somewhat softer episode (lasting ~100 s). 
The total energy content of the second episode is 5 times that of the first. Its 
light curve shape is similar in many respects to that of short bursts with ex- 
tended emission. There are, however, differences in that the short episode of this 
event is longer than the previously detected examples and the soft episode is rel- 
atively brighter. Another similarity with short bursts comes from a lag analysis 
of GRB 060614. The lag between temporal structures in the 50 — 100 keV band 
and those in the 15 — 25 keV bands for the first 5 s is 3 it 6 ms which falls in the 
same region of the lag-luminosity plot as short bursts (Figure 5). It is difficult 
to determine unambiguously which category of burst GRB 060614 falls into. It 
is a long event by the traditional definition, but it lacks an associated SN as 
had been observed in all other nearby long GRBs. It shares some similarities 
with Swift short bursts, but has important differences such the brightness of the 
extended soft episode. If it is due to a collapsar, it is the first indication that 



some massive star collapses either fail as supernovae (Woosley 19931 or highly 



underproduce ^^Ni ( Lopez-Camara et al.||2008 |; if it is due to a merger, then the 
bright long-lived soft episode is hard to explain within the framework of compact 



binary mergers (Di Matteo et al. 2002 Lee et al. 2004, Setiawan et al. 2004 1 



Thus, this peculiar burst has challenged the usual classifications of GRBs. 



5 BASIC PHYSICAL CONSIDERATIONS 

In this section, we endeavor to outline some of the physical processes that are 
believed to be most relevant to interpreting GRBs. Though the field is far from 
maturity, sufficient progress has been made in identifying the essential physical 
ingredients. A basic scheme can provide a conceptual framework for describing 
the observations even when the framework is inaccurate! The following should 
be interpreted in this spirit. 
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Figure 14 Diagram exhibiting GRB activity over successive decades in radius ranging 
from 10^ cm to Ipc. The first task in attempting to construct a general scheme of 
GRBs is to decide which parameters exert a controlling influence upon their properties. 
(10^ cm) The GRB nucleus (BH or NS) and its magnetosphere. (10^ cm) The accretion 
flow is likely to be embedded in a very active corona. We expect coronal arches, as well 
as large magnetic structures, to be quite common and to be regenerated on an orbital 
timescale. Relativistic outflow from the black hole is proposed to be focused into two 
jets. (10^*^ cm) Even if the outflow is not narrowly coUimated, some beaming is expected 
because energy would be channelled preferentially along the rotation axis. The majority 
of stellar progenitors, with the exception of some very compact stars, will not collapse 
entirely during the typical duration of a GRB. A stellar envelope will thus remain to 
impede the advance of the jet. (10^^ cm) This is the typical size of an evolved massive 
star progenitor. A thermal break-out signal should precede the canonical, softer 7-rays 
observed in GRBs. (10^^ cm) Velocity differences across the jet proflle provide a source of 
free energy. The most favourable region for shocks producing highly variable 7-ray light 
curves is above the baryonic or pair-dominated photosphere. (10^^ cm) The external 
shock becomes important when the inertia of the swept-up external matter starts to 
produce an appreciable slowing down of the ejecta. (10^^ cm) Finally, we come to the 
end of the relativistic phase. This happens when the mass E /c^ has been swept-up. 



GRB activity manifests itself over a dynamical range of ~ 13 decades in ra- 



dius. In Figure 14 we show a schematic montage of successive decades, exhibiting 
phenomena which are believed to take place on each of these length scales. The 
phenomena are not directly observed and the associated frames represent edu- 
cated guesses of their geometrical arrangements. An anatomical summary of the 
underlying physical processes working outward from the smallest to the largest 
scales follows. 

5.1 The Central Engine 



In principle, flow onto a compact object can liberate gravitational potential en- 
ergy at a rate approaching a few tenths of Mc^, where M is the mass inflow 
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rate. Even for such high efficiencies the mass requirements are rather large, with 
the more powerful GRB sources ~ 10^'^ erg s~^ having to process upwards of 
IO^^Mq s^"*^ through a region which is not much larger than the size of a NS 
or a stellar mass BH. Radiation of the BH rest mass on a timescale rg/c, where 
Tg = GM/(? = 1.5 X 1O^(M/M0) cm is the characteristic size of the collapsed 
object, would yield luminosities c' /G = 4 x 10^^ erg s~^. 

When mass accretes onto a BH or NS under these conditions, the densities and 
temperatures are so large that photons are completely trapped and neutrinos, 
being copiously emitted, are the main source of cooling. The associated interac- 
tion cross section is then many orders of magnitude smaller, and, as a result, the 
allowed accretion rates and luminosities are correspondingly much higher. For 
example, using the cross section for neutrino pair productioij^ the Eddington 
limit can be rewritten as LEdd,;. = 8 x lO^^(£^,,/5OMeV)"^(M/M0) erg s'S with 
an associated accretion rate is MEdd,i/ x (efficiency)"^, where MEdd,u = -^^Edd,i//c 
( Ramirez-RuizpOOGb Lee and Ramirez-Ruiz|[2007| ). 



The blackbody temperature if a luminosity L-Edd,v emerges from a sphere of 
radiu^Tg is 

The radiation temperature is expected to be < Tth = GMmp/(3krg) ~ 200 MeV, 
the temperature the accreted material would reach if its gravitational potential 
energy were turned entirely into thermal energy. 

Related to this, there is a fiducial density in the vicinity of the BH 

= ^ ~ 10"(M/Me)-' , cn,-. (2) 

It should be noted that the typical Thomson optical depth under these conditions 
is tt ~ n^^d u'^z ~ ^^^^ expected, photons are incapable of escaping 

and constitute part of the fiuid. A characteristic magnetic field strength is that 
for which B^/Svr = riEdd,!/ "ipC^: 

- ("w) - ' ^ ^iMiM^^n (_^) a (3) 

If a field -BEdd,!/ were applied to a BH with J « Jmax, the electromagnetic power 
extraction would be ~ L-p^ddM (Usov|[T992 1. 



5.2 Accretion Flows 

As we discussed above, a BH or NS embedded in infalling matter offers a more 
efficient power source than any other conceivable progenitor. Although this ef- 
ficiency is over a hundred times larger than that traditionally associated with 

^Although we have considered here the specific case of neutrino pair creation, the estimates 
vary little when one considers, for example, coherent scattering of neutrinos by nuclei and/or 
free nucleons (except for the energy scaling). 

^Similar overall fiducial numbers also hold for neutron stars, except that the simple mass 
scalings obtained here are lost. 
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thermonuclear reactions, the required rate of mass supply for a typical GRB is 
of course colossal. Such high mass fueling rates are never reached for BHs in 
XRBs or AGN, where the luminosity remains well below the photon Eddington 
limit. They can, however, be achieved in the process of forming NS and stellar- 



mass BHs during the collapse of massive stellar cores (Houck and Chevalier||1991 


Woosley||1993 MacFadyen and Woosley|1999 Narayan et al.|2001 


) and in binary 


mergers involving compact objects (Ruffert et al. 


1997 Kluzniak and Lee 1998 


Rosswog and Liebendorfer 2003 , Lee et al. 2005b 


Setiawan et al. 


2006 Metzger 



been directed towards describing the possible formation channels for these sys- 



tems (Fryer et al. 1999, Belczynski et al. 2002, Izzard et al. 2004 Podsiadlowski 



et al. 



2004). This involves evaluating those which are likely to produce a vi- 



able central engine, and understanding the flow patterns near relativistic objects 



accreting matter in the hypercritical regime (Di Matteo et al. 2002 Chen and 



Beloborodovl|2007| |Proga and Zh^|2006l |Perna et al.|[2006l [Shibata et al]|2007 



Kumar et al. 20081, where photons are unable to provide cooling but neutrinos 
do so efficiently. 

The angular momentum is quite generally a crucial parameter, in many ways 



determining the geometry of the flow (Lee and Ramirez- Ruiz 20061. The quasi- 
spherical approximation breaks down when the gas reaches a radius rcirc ~ 
I'^/GM, where I is the angular momentum per unit mass, and if injection oc- 
curs more or less isotropically at large radii, a familiar accretion disk will form. 
The inner regions of disks with mass fluxes < MEdd,;/ generally able to cool 
by neutrinos on time scales shorter than the inflow time. 

If m = M/M^dd,u ^ 1) then the bulk of the neutrino radiation comes from 
a region only a few gravitational radii in size, and the physical conditions can 
be scaled in terms of the Eddington quantities deflned above. The remaining 
relevant parameter, related to the angular momentum, is finflow/f^freefaib where 
Wfreefaii — (2GM/i?)^/^ is the free fall velocity. The inward drift speed i^infiow would 
be of order Ufreefall for supersonic radial accretion. When angular momentum is 
important, this ratio depends on the mechanism for its transport through the 
disk, which is related to the effective shear viscosity. For a thin disk, the factor 
(finflow/'yfreefaii) is of Order a{H / R)'^ , where H is the scale height at radius R and 



a is the phenomenological viscosity parameter (Shakura and Syunyaev 19731. 

The characteristic density, at a distance r from the hole, with account of the 
effects of rotation, is 



p ~ m(r/rg) ^/^(Wmflow/^^freefalOPEdd,!., 



(4) 



and the maximum magnetic fleld, corresponding to equipartition with the bulk 
kinetic energy, would be 



Bcq ~ m^/^(r/rg) ^/"^(I'mflow/'WfrcefalO^^^-B: 



Edd,i/- 



(5) 



Any neutrinos emerging directly from the central core would have energies of a 
few MeV. Note that, as mentioned above, kTEdd,u is far below the virial tem- 
perature /cTvir — mpC^{r /rg). The flow pattern when accretion occurs would be 
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then determined by the value of the parameters Liy/LEdd,u, which determine the 
importance of radiation pressure and gravity, and the ratio tcooi/^dynamicai) which 
fixes the temperature if a stationary flow pattern is set up. 

It is widely believed that accretion onto a compact object, be it a NS or a stellar 
mass BH, offers the best hope of understanding the "prime mover" in all types of 
GRB sources although a possible attractive exception includes a rapidly spinning 



neutron star with a powerful magnetic field ( 


Usov 


1992 


Thompson||1994 


Spruit 


1999 [Thompson et al. |2004 Uzdensky and MacFadyen 2007 Bucciantini et al. 



2008 



Dessart et al. 2008 1 . If there is an ordered field B, and a characteristic 
angular velocity lo, for a spinning compact source of radius r^,, then the magnetic 



dipole moment is ~ Br'^. General arguments suggest (Gold 1968 Pacini 19671 
that the non-thermal magnetic-dipole-like luminosity will be Lem ~ -B^r^w^/c^, 
and simple scaling from these familiar results of pulsar theory require fields of 
order 10^^ G to carry away the rotational or gravitational energy (which is ~ 
-|^q53 gj.g^ jj-^ time scale of seconds ( Usov|[l992 1. 



5.3 Jet Production 

There are two ingredients necessary for the production of jets. First, there must 
be a source of material with sufficient free energy to escape the gravitational 
field of the compact object. Second, there must be a way of imparting some 
directionality to the escaping flow. Our eventual aim must be to understand 
the overall flow pattern around a central compact object, involving accretion, 
rotation, and directional outflow but we are still far from achieving this. Most 
current workers who have discussed outflow and coUimation have simply invoked 
some central supply of energy and material. A self-consistent model incorporating 
outflow and inflow must explain why some fraction of the matter can acquire a 
disproportionate share of energy (i.e., a high enthalpy). 

A spinning black hole (or neutron star) constitutes an excellent gyroscope, 
and the ingredients of accretion, angular momentum, entropy production, and 
possibly magnetic fields are probably sufficient to ensure the production of coUi- 



mated outflow (McKinney 2006). However, the detailed mechanism is a matter 
for debate (it is not even clear what is accelerated), and several distinct lines 
of research are currently being pursued. One solution is to reconvert some of 
this energy via collisions outside the disk into electron-positron pairs or photons 



(Mochkovitch et al.|1993 


Goodman et al.|1987 


Rosswog and Ramirez-Ruiz|2002 


Rosswog et al. 


2003 


Aloy et al. 


2005 


Dessart et al. 


2009 


). If this occurs in a 



region of low baryon density (e.g., along the rotation axis, away from the equa- 
torial plane of the disk) a relativistic pair-dominated wind can be produced. An 
obvious requirement for this mechanism to be efficient is that the neutrinos es- 
cape (free streaming, or diffusing out if the density is high enough) in a time scale 
shorter than that for advection into the BH. The efficiency for conversion into 
pairs (scaling with the square of the neutrino density) is too low if the neutrino 
production is too gradual, so this can become a delicate balancing act. 

Jets may alternatively be produced electromagnetically. Such mechanism could, 
in principle, circumvent the above restriction in efficiency. The potential differ- 
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ence across a disk threaded by an open magnetic field lines can exceed 10^2 y, 
and this is available for accelerating high-energy particles which will produce 
an electron-positron cascade and ultimately a relativistic jet that carries away 
the binding energy of the accreting gas. Blandford & Znajek (1977) extended 
this idea to black holes and showed how the spin energy of the hole could like- 
wise be extracted. A hydromagnetic description of this mechanism is more likely 



to be appropriate (Blandford and Payne 1982 1 . The field required to produce 
Lem ^ 10^^ erg s~^ is enormous, and may be provided by a helical dynamo op- 



erating in hot, convective nuclear matter with a millisecond period (Duncan and 



Thompson 1992 1. A dipole field of the order of 10^^ G appears weak compared 



to the strongest field (Price and Rosswog 2006 Thompson et al. 20051 that can 



in principle be generated by differential rotation (~ 10^^[P/1 ms]~^ G), or by 
convection {^^ 10"*^^ G), although how this may come about in detail is not re- 
solved. Note, however, that it only takes a residual torus (or even a cold disk) of 
10-3 Mq to confine a field of 10^^ G. 

A potential death-trap for such relativistic outflows is the amount of entrained 
baryonic mass from the surrounding medium. For instance, a Poynting flux of 
■j^q53 gj.g (.Q^if^j ]2ot accelerate an outflow to P > 10^ if it had to drag more than 
~ IO-^Mq of baryons with it. A related complication renders the production 
of relativistic jets even more challenging, because the high neutrino fluxes are 
capable of ablating baryonic material from the surface of the disk. Thus a rest 
mass flux limits the bulk Lorentz factor of the wind to P,; = L^^i^^/ [MrjC^). 
Assuming that the external poloidal field strength is limited by the vigor of 
the convective motions, the spin-down luminosity scales with neutrino flux as 
^wind ~ -^^em Gc cc v'^ ^^ PC L^^^, where Vcon is the convective veloci ty. The 
^/^ /iiv/Tof rjrroT- of oi onnv o^on \KTr^r^^^^^r ioo«l\ whlch 



ablation rate is oc Lj (Metzger et al. 



2007 



Qian and Woosley 



19961 



indicates that the limiting bulk Lorentz factor P^ of the wind decreases as L~^. 
Thus the burst luminosity emitted by a magnetized neutrino cooled disk may be 
self-limiting. Mass loss could, however, be suppressed if the relativistic wind were 
somehow collimated into a jet. This suggests that centrifugally-driven mass loss 
will be heaviest in the outer parts of the disk, and that a detectable burst may 
be emitted only within a relatively small solid angle centered on the rotation axis 
( Levinson and Eichler||2000 ) . 



5.4 Jet Collimation, Stability, and Confinement 

As we discussed above, one of the key issues one needs to address concerning jets 
is how they are formed? A second one is, how do they retain their coherence 
and collimation as they transverse circumburst space? The second issue is no 
less pressing than the first, since jet-like flows known on Earth are notoriously 
unstable. The situation in a jet is nonetheless different from that encountered in 
the laboratory because of the super- Alfvenic and supersonic streaming velocity 



(McKinney and Blandford 20081. MHD probably provides a better description 



of the macroscopic behavior of a GRB jet than in the case of laboratory plasmas 
because the particle Larmor radii are so much smaller than the transverse size 
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of the jet. Formal stability analysis of even the simplest jet models is complex. 
Unstable MHD modes (such as pinch and kink instability) do not grow as rapidly 
as in stationary plasma, although sufficiently short-wavelength unstable modes, 
localized within the jet interior are not suppressed by the relative motion. Lon- 
gitudinal magnetic field in the core of the jet can likewise act as a backbone, 
provided that the correlation length of the magnetic field reversal along the jet 
exceeds the wavelength of the perturbation. Instabilities may be suppressed for 
perturbations of large wavelength along the jet by the inertia of the ambient 
medium. Whatever one's view of the relative merits of fiuid and electromagnetic 
models of relativistic jets - and perhaps the truth lies between the two extremes 
- it is clear that our understanding of the role of magnetic fields in jets is less 
advanced than that of other aspects of the problem. 

An understanding of the collimation and confinement of a jet can come about 
only through a knowledge of the properties of the medium through which it 
propagates. Information about the ambient pressure gradient propagates into 
the jet at the internal sound speed. If the jet moves through a background with 
pressure scale height ~ r, the necessary condition for the jet interior to remain in 
pressure balance with its surroundings is 9Mj < 2, where Mj is the Mach number 
and ~ d/r is the opening angle. If this condition is not satisfied when a jet is 
moving into a region of higher pressure, then strong shocks are driven into the 
jet. A jet moving into a region of lower pressure becomes overpressured relative 
to its surroundings, and thereafter expands freely. 

Freedom or confinement? This is the first question for which an extrinsic or 
environmental effect comes into play (which may, in turn, strongly affect what 
we observe), and it is thought to be particularly important for massive stellar 
GRB progenitors. This is because in such stars, a stellar envelope will remain 



to impede the advance of the jet ( 


MacFadyen and Woosley 


1999 


MacFadyen 


et al.||2001 Meszaros and Rees||2001 Aloy et al.||2000 Ramirez-Ruiz et al.||2002a| 


Matzner 


2003 


Zhang et al. 


2004 


Proga et al. 


20031. The beam will then evac- 



uated a channel out to some location where it impinges on the stellar envelope. 
A continuous flow of relativistic fluid emanating from the nucleus supplies this 
region with mass, momentum, energy, and magnetic flux. Most of the energy 
output during that period is deposited into a cocoon or wastehasket surround- 
ing the jet, which, after expansion, would have enough kinetic energy content to 
substantially alter the structure of the relativistic outflow, if not, in fact, provide 
much of the observed explosive power (Ramirez-Ruiz et al. 2002a I . 

Can relativistic jets really be formed inside stars? This is not the most pro- 
pitious environment to create an ultra relativistic, baryons-starved jet. What 
needs to be demonstrated, is that the outflow is not "poisoned" by baryons by 
the time that it reaches the surface of the star. It appears that it is not necessary 
to coUimate the jet very tightly or to achieve high bulk Lorentz factor as the 
flow leaves the stellar surface (Zhang et al. 2004 Morsony et al. 2007). As long 
as the emergent flow has a high enthalpy per baryon, it will expand and achieve 



its high terminal speed some distance from the star (Uzdensky and MacFadyen 



2007 Bucciantini et al. 20081. A strong thermal break-out signal is expected to 
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precede the canonical 7-rays observed in GRBs with massive progenitors as the 
shock breaks through the stellar surface and exposes the hot shocked material 



Waxman et al. 



(MacFadyen et al. 2001 Ramirez- Ruiz et al. 2002c, Waxman and Meszaros 2003 



2007). For very extended envelopes, the jet may be unable to 
break through the envelope. TeV neutrino signals produced by Fermi acceler- 
ated relativistic protons within the cork may provide a means of detecting such 
choked-off, 7-ray dark collapses ( Meszaros and Waxman||2001 1. 



5.5 Dissipation and Cooling Effects Within the Jet 

The unique feature of GRBs is that the bulk Lorentz factor T may reach values 
from hundreds to thousands ( Lithwick and Saji]|2001 1. The relativistic motion of 
the radiating particles introduces many interesting effects in GRB emissions that 
must be properly taken into account (e.g. Fenimore et al.|[l996b |. 

Three frames of reference are considered when discussing the emission from 
systems moving with relativistic speeds: the stationary frame, which is denoted 
here by asterisks, the comoving frame, denoted by primes, and the observer frame. 
The differential distance traveled by the expanding source during differential time 
dt* is dr = Pcdt^, where /3 = Vl — r~2. Due to time dilation, dr = (JTcdt' . The 
relationship between comoving and observer times is (1 -|- z)Tdt'{l — (3 cos 6) = 
(1 -|- z)dt' /6 = dt, where 9 is the angle between the emitting element and the 
observer, 6 = [r(l — /3cos^)]~^ is the Doppler factor and z the cosmological 
redshift. For an on-axis observer we therefore see that dt = {1 + z)dr /T'^c, and, as 
a result, the blast wave can travel a large distance T^cAt during a small observing 



hv/meC^ 



IS 



time interval. A photon detected with dimensionless energy e 
emitted with energy 5e' /{I + z). 

Few would dispute the statement that the photons which bring us all our 
information about the nature of GRBs are the result of triboluminescence. For 
instance, velocity differences across the jet profile provide a source of free energy 
from particle acceleration through shock waves, hydromagnetic turbulence, and 



tearing mode magnetic reconnection (Rees and Meszaros 1994 1 . If the value of 



r at the base increases by a factor > 2 over a timescale At, then the later ejecta 
will catch up and dissipate a fraction of their energy at radius given by 



cAtr^ ~ 3 X lO^^AtoF^ cm 



(6) 



where At = lAto s and F = 10 Dissipation, to be most effective, must occur 
when the wind is optically thin: tt — n'(JT{r/T) < 1 (here n' is the comov- 
ing number density). Otherwise it will suffer adiabatic cooling before escaping 



(Goodman 1986). The photosphere (baryonic or pair-dominated) is a source of 



soft thermalized radiation, which may be observationally detectable in some GRB 
spectra and may also result in inverse Compton cooling of the nonthermal elec- 



trons accelerated in the shocks occurring outside it (Meszaros and Rees 2000 



Ramirez-Ruiz and Lloyd-Ronning 2002 , Sp ada et al.||2000[ Kobayashi et al.||2002 



Meszaros et al. 2002', 'Ram irez-RuiZ||2005t |Giannio s and Spruitj 2007 Pe'er et al, 



2006[ |Ryde and Pe'er||2008( [Thompson et al.„2007| ). 
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In the presence of turbulent magnetic fields built up behind the shocks, the 
electrons can produce a synchrotron power-law radiation spectrum, whereas the 
inverse Compton scattering of these synchrotron photons extends the spectrum 



into the GeV range (Meszaros et al. 1994 1 . To illustrate the basic idea, suppose 



that electrons, protons, and magnetic field share the available internal energy, 
then the electrons reach typical random Lorentz factors of 7 ~ mp/me] while 
the assumption of a Poynting flux L b implies a comoving magnetic field of order 



B 



rl/2 
ll3 



-1 ~ lO^L^/gor^^i^ra^ G, where Lb = 10^^Lb,5o erg s^^ and 
W^'^r^^i^. For these values of 7 and B, the typical observed synchrotron 
frequency is Usy ~ 0.5L^J'^^Qr~l^{l + z)~^ MeV, independent of the bulk Lorentz 
factor r, and in excellent agreement with the observed values of the uF^^ peak 
of GRB spectra (Figure 1). Yet, there are in some instances serious problems 
associated with this model (e.g., dissipation efficiencies). These difficulties have 
motivated consideration for alternative scenarios (e.g., Kumar and Narayan|2008 
Kumar and McMahon||2008|). 



A magnetic field can ensure efficient cooling even if it is not strong enough 
to be dynamically significant. If, however, the field is dynamically significant 



in the wind (Meszaros and Rees 19971, then its internal motions could lead to 
dissipation even in a constant velocity wind (Thompson 1994 1 . Instabilities in 



this magnetized wind may be responsible for particle acceleration (Thompson 



2006), and it is possible that 7— ray production occurs mainly at large distances 



from the source (Lyutikov and Blandford 2003) 



A further effect renders the task of simulating unsteady winds even more chal- 
lenging. This stems from the likelihood that any entrained matter would be a 
mixture of protons and neutrons. If a streaming velocity builds up between ions 
and neutrons then interactions can lead to dissipation even in a steady jet where 
there are no shocks dPerishev et"aL]|1999[ [B'eloborodov||2003[ ). 



5.6 Jet Interaction with the External Environment 

Astrophysicists understand supernova remnants reasonably well, despite contin- 
uing uncertainty about the initiating explosion; likewise, we may hope to under- 
stand the afterglows of GRBs, despite the uncertainties about the trigger that 
we have already emphasized. 

In the simplest version of the afterglow model, the blast wave is approximated 
by a uniform thin shell. A forward shock is formed when the expanding shell 
accelerates the external medium, and a reverse shock is formed due to deceleration 
of the cold shell. The forward and reverse shocked fiuids are separated by a 
contact discontinuity and have equal kinetic energy densities. 

As the blast wave expands, it sweeps up material from the surrounding medium 



to form an external shock (Meszaros and Rees 1993). Protons captured by the 



expanding blast wave from the external medium will have total energy TnipC 
in the fiuid frame, where rup is the proton mass. The kinetic energy swept into 
the comoving frame by an uncoUimated blast wave at the forward shock per unit 



time is given by ( [Blandford and McKee||l976[ ) dE'/dt' = 4TTr^neKtmpC^pr{r - 1) 
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where factor of T represents the increase of external medium density due to length 
contraction, the factor (F — 1) is proportional to the kinetic energy of the swept 
up particles, and the factor P is proportional to the rate at which the particles 
are swept. 

The external shock becomes important when the inertia of the swept-up ex- 
ternal matter starts to produce an appreciable slowing down of the ejecta. The ex- 
panding shell will therefore begin to decelerate when E = TMbC^ = r^mpC^(47rrjjn 
giving the deceleration radius ( Meszaros and Rees|1993 Rees and Meszaros|1992 1 



3E 



1/3 



rd 



3 X 10 



16 



Vr^ncxt/ 



1/3 



cm, 



(7) 



where T = E/My^c^ is the coasting Lorentz factor, Mb is the baryonic mass, 
E = 10^^ £"52 erg is the apparent isotropic energy release and next is the number 
density of the circumburst medium. This sets a characteristic deceleration length. 
This deceleration allows slower ejecta to catch up, replenishing and re-energizing 
the reverse shock and boosting the momentum in the blast wave. 

Most treatments employing blast-wave theory to explain the observed after- 
glow emission from GRBs assume that the radiating particles are electrons. The 
problem here is that ~ nip /me ~ 2000 of the nonthermal particle energy swept 
into the blast-wave shock is in the form of protons or ions, unless the surround- 
ings are composed primarily of electron-positron pairs. For a radiatively efficient 
system, physical processes must therefore transfer a large fraction of the swept- 



up energy to the electron component (Gedalin et al. 2008). In most elementary 



treatments it is simply assumed that a fraction eg of the forward-shock power is 
transferred to the electrons. 

The strength of the magnetic field is another major uncertainty. The stan- 
dard prescription is to assume that the magnetic field energy density is a 
fixed fraction of the downstream energy density of the shocked fluid. Hence 
i?V(8vr) 



UB 



r) (although see, e.g. 



Rossi and Rees 



2003 ) . It is also generally supposed in simple blast- wave model calculations that 



some mechanism injects electrons with a power-law distribution between elec- 
tron Lorentz factors 7min < 7 < 7max downstream of the shock front, where the 
maximum injection energy is obtained by balancing synchrotron losses and an 
acceleration rate given in terms of the inverse of the Larmor time scale. 

A break is formed in the electron spectrum at cooling electron Lorentz factor 
7c, which is found by balancing the synchrotron loss time scale tgy with the adi- 
abatic expansion time t^^j 
wave, r oc t~^/^, so that 7n 



r/{Tc) (Sari et al. 19981. For an adiabatic blast 
oc t~^/^ and 7c oc t^l^. As a consequence, the ac- 
celerated electron minimum random Lorentz factor and the turbulent magnetic 



field also decrease as inverse power-laws in time (Meszaros and Rees 1993 Rees 



and Meszaros|[T992 1 . This implies that the spectrum softens in time, as the syn- 
chrotron peak corresponding to the minimum Lorentz factor and field decreases. 



leading to the possibility of late long wavelength emission (Sari et al. 1998). 



The relativistic expansion is then gradually slowed down, and the blast wave 
evolves in a self-similar manner with a power-law lightcurve. This phase ends 
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when so much mass shares the energy that the Lorentz factor drops to 1 (Ayal 



and Piran 2001, Ramirez-Ruiz and MacFadyen 2008). Obviously, this happens 



when a mass E/c has been swept up. This sets a non-relativistic mass scale: 



p2/3 



3^ 



47rmpc2ncxt 



\ 1/3 

2. 1 



10 



18 



next/ 



cm 



(8) 



For comparison, the Sedov radius of a supernova that ejects a 10 Mq envelope 
could reach ~ 5 pc or more. 

In GRB sources, with jets which we believe to be highly relativistic, the orienta- 
tion of the jet axis with respect to our line of sight will strongly affect the source's 
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|l995l 
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Granot 
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2007), because radiation from jet material will be Doppler 
beamed in the direction of motion. Attempts to understand the luminosity func- 
tion of GRBs may have to take into account the statistics of orientation, collima- 
tion,and velocity of the jet, as well as the jet's intrinsic radiation properties. 



Although our proposed synthesis of GRB physical properties is highly conjec- 
tural and far from unique, we hope it will provide a framework for discussing 
the integrated properties of these objects. We conclude the review by discussing 
how future observations, experiments and theoretical studies should enhance our 
understanding of the physical properties underlying GRBs. 



6 PROSPECTS 



GRB studies, especially the afterglow-enabled studies of the last ten years, re- 
main a young field. The years ahead are thus sure to bring astonishing discoveries 
as the capabilities and experience of observers improve, theorists make more and 
stronger ties to physical theory, and new and upgraded facilities open vistas. In 
this section, we summarize the instrumental capabilities and theoretical oppor- 
tunities for near-term progress in GRB research. 



6.1 Facilities 



Table [2] summarizes recent and near-term future GRB missions, planned and 



proposed. Swift (Gehrels et al. 2004) is now the primary mission and has excellent 



prospects for continued operation, with an orbit that will be stable until at least 
2020. Fermi and AGILE, also in continuing operations, are providing added burst 
detections with simultaneous high energy (>100MeV) coverage that promises to 
redefine the maximum energies and Lorentz factors that GRB engines are capable 
of producing. 

The SVOM mission ( Paul et al.||2008" ), currently under development, promises 
4 to 300 keV coverage and Swift-like slews that will bring X-ray and optical tele- 



scopes to bear on burst positions. The proposed JANUS small explorer (Roming 



2008) would focus on high-redshift bursts, detecting prompt emission over the 1 



to 20 keV band and slewing to observe afterglows with a near-infrared telescope 
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Table 2 Recent and future GRB Missions 



Mission 


Trigger energy range 


FOV 


Detector area 


Other wavelengths 


GRB rate (yr ) 


BATSE 


20keV-1.9 MeV (LAD) 
10 keV-100 MeV (oD) 


4tt sr 


2025 cm^ per LAD 
127 cm per bD 




300 


HETE-2 


6-400 keV 


3 sr 


120 cm^ 


X-ray 




Swift 


15-150 keV 


1.4 sr 


5200 cm2 


UV, Opt, X-ray 


100 (~10% SGRBs) 


AGILE 


30 McV-50 GeV 


~3 sr 




Hard X-ray 




Fermi 


20 MeV-300 GeV 
8 keV-1 MeV 
150 keV-30 MeV 


>2 sr (LAT) 
9.5 sr (GBM) 


>8000 cm2 (LAT) 
126 cm2 (GBM-LED) 
126 cm2 (GBM-HED) 




50 


SVOM 


4 keV-300 keV 
50 keV-5 MeV 


2 sr (CXG) 

89° X 89° (GRM) 




Optical, X-ray 


80 


JANUS 


1-20 keV 


4 sr 




Near-IR 


25 (high z) 


EXIST 


5-600 keV 


~3.6 sr (HET) 


5.96 m2 (HET) 


Optical, near-IR, 
X-ray 


300 



Note: LAD, Large Area Detector; SD, Spectroscopy Detector; LAT, Large Area Telescope; LED, Low Energy Detector; HED, High 
Energy Detector; CSG, X-ray/Gamma-ray Camera; GRM, Gamma-Ray Monitor; GBM, GLAST Burst Monitor; SGRB, short GRB. 



(50 cm aperture, 0.7 to 1.7 /um coverage). The proposed large mission EXIST has 
also added fast-response slews and focusing X-ray and near-infrared telescopes to 
its original complement of hard X-ray imaging detectors ( |Grindlay|[2008| ) . 

Current and near-future ground-based facilities include a fast-growing array of 
robotic telescopes primed to respond to burst alerts, ongoing improvements to the 
instrumentation and capabilities of large-aperture telescopes, a new generation 
of air-Cerenkov TeV facilities, and with the EVLA initiative, the most significant 
upgrade to a high- sensitivity radio facility in decades. 

Rapid follow-up of GRB discoveries is occurring both in space and on the 
ground; the Swift UVOT observes ~80 GRBs per year within 2 minutes, detect- 
ing about half of them. On the ground, fast new telescopes such as ROTSE-III 
dAkerlof et ar][2003| ), RAPTOR ( [Wozniak et al][2006l ), Pi-of-the-sky ( |Burd et al 



2005D , and REM ( |Zerbi et al.||2001[ ) are able to point at GRBs within 10 - 20 s. 



and larger facilities like the seven-band GROND ( Greiner et al.|2008 ) respond on 
10-minute timescales. Ready availability of sensitive CCDs and HgCdTe detec- 
tors, combined with rapid-slew mounts and autonomous software systems, should 
enable further expansion of rapid-response telescopes with more larger (Z) ^ 2 m) 
facilities anticipated in the near future. 

GRB-related programs continue to compete successfully for time on premier 
optical facilities. Large-aperture (D ;^ 6 m) telescopes provide the spectroscopic 
observations necessary for GRB redshift measurements. At the VLT, a rapid 
response mode has gone further to provide time-sequence observations revealing 
variable absorption from the host galaxy, and the impending commissioning of 
the X-Shooter spectrograph ( Kaper et al. 2009 1 will enable full UV to NIR charac- 
terization of afterglow spectra with a single integration. In addition to the TOO 
opportunities provided by GRB alerts, multiple host galaxy survey programs are 
underway. 



At radio wavelengths, the workhorse facilities have been the VLA (e.g., Chan- 



dra et al. 20081 and WSRT (e.g., van der Horst et al. 2007), providing data 



primarily in the 1 — 10 GHz at sensitivities in the 0.2 — 1.0 mJy range; the 
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Expanded VLA (EVLA) upgrade will improve the sensitivity of that facility to 
'^lO /iJy. Looking ahead, by 2012, the ALMA array will be on line operating in 
the higher frequency range 90 — 950 GHz with > 100 times the sensitivity of the 
VLA. The peak in the synchrotron spectrum for a wide range of GRB parame- 
ters lies in the ALMA range, and should make it a powerful future tool for radio 
observations. 

Observations of GRBs at TeV energies can be performed by both narrow- 
field air Cerenkov facilities and wide- field water Cerenkov detectors; the latter 
approach, at Milagro, yielded the tentative detection of prompt TeV emission 



for GRB 970417A (Atkins et al. 2003 1. There has been no detection to-date 



of prompt or afterglow emission with the air Cerenkov facilities, but significant 
effort in pursuit of burst alerts is underway at MAGIC, HESS, and VERITAS. 
The MAGIC dish, in particular, has a rapid-response mode that has provided 



observations of multiple bursts within a minute of trigger ( Albert et al. 2006 1 . 
Future facilties will seek to lower detection thresholds to ~ 100 GeV, which 
would provide a significantly-expanded horizon within which GRB sources will 
be visible, rather than attenuated by photon-photon interactions. 

In space, the upcoming HST SM4 promises to revive that facility, not only 
providing new and resuscitated instruments, but also the gyros necessary for 
flexible and fast-response scheduling. In the future, the capabilities of JWST 
will provide an excellent resource for high signal-to-noise observations of highly 
obscured and high-redshift GRBs. 

Within the X-ray band, the Swift XRT has redefined all expectations for the 
characterization of X-ray afterglows, and at the same time, proved a highly ef- 
fective facility for refining the multi-arcmin localizations provided by other GRB 
missions. Chandra, XMM, and Suzaku continue to be active in GRB observations, 
and future missions including the Indian ASTROSAT and Japan/US Astro-H 
promise added capabilities in the near future. Dramatic improvements will await 
the next-generation GRB facilities (Table |2]) or the arrival of the International 
X-ray Observatory. 

6.2 Multimessenger Aspects 

Given the rapid ongoing expansion of the capabilities of non-electromagnetic de- 
tector facilities, and the extreme luminosity and time- specificity of GRB sources, 
it will not be long before we either have the first coincident, multi- messenger de- 
tection of a GRB, or realize limits on the non-electromagnetic emissions of GRBs 
that challenge our current understanding. 

The same shocks which are thought to accelerate electrons responsible for the 
non-thermal 7-rays in GRBs should also accelerate protons, leading ultimately to 
copious emission of high-energy neutrinos ( Waxman|2004b Waxman|2004c Wax- 



man ] [2006l [Dernier and Holmes]|2005l ) . The maximum proton energies achievable 



in GRB shocks are Ep ~ 10^*^ eV, comparable to the highest energies measured 
with large cosmic ray ground arrays ( Hayashida et alT]|1999 |. For this, the accel- 



eration time must be shorter than both the radiation or adiabatic loss time and 



the escape time from the acceleration region (Waxman 19951. The accelerated 
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protons can interact with the fireball photons, leading to charged pions, muons 
and neutrinos. For internal shocks producing observed 1 MeV photons this im- 
plies > 10^^ eV protons, and neutrinos with ~ 5% of that energy, > 10^^ eV 
(Waxman and Bahcall 1997). Another copious source of target photons in the 
UV is the afterglow reverse shock, for which the resonance condition requires 
higher energy protons leading to neutrinos of 10^^ — 10^^ eV (Waxman and Bah- 



call 1999). Whereas photon-pion interactions lead to higher energy neutrinos 



and provide a direct probe of the shock proton acceleration as well as of the pho- 
ton density, inelastic proton-neutron collisions may occur even in the absence of 
shocks, leading to charged pions and neutrinos ( Derishev et al.|[l999 ) with lower 
energies than those from photon-pion interactions. The typical neutrino ener- 
gies are in the ^ 1-10 GeV range, which could be detectable in coincidence with 
observed GRBs. This is the province of projects like AMANDA, IceCube and 
ANTARES. Neutrino astronomy has the advantage that we can see the universe 
up to ~ EeV energies. By contrast, the universe becomes opaque to 7-rays above 
~ TeV energies through absorption by the infrared background. 

The last and most challenging frontier is that of gravitational radiation, which 
is largely unknown territory. A time-integrated luminosity of the order of a solar 
rest mass (~ 10^^ erg) is predicted from progenitor models involving merging 
compact objects, while that from coUapsar models is less certain, and expected 
to be lower by at least one order of magnitude. 

Ground-based facilities, like LIGO, TAMA and VIRGO, are currently seeking 
the first detection of these stellar-scale, high-frequency {v ^ 50 Hz) sources. The 
observation of associated gravitational waves would be facilitated if the merg- 
ers involve observed short GRB sources; and conversely, it may be possible to 
strengthen the case for (or against) NS-NS or NS-BH progenitors of short bursts 
if gravitational waves were detected (or not) in coincidence with some events. 
The technical challenge of achieving the sensitivities necessary to measure waves 
from assured sources should not be understated; neither, however, should the 
potential rewards. 

The Enhanced LIGO interferometers will be on line in 2009 with ability to 
detect NS binary mergers to 20 Mpc. The Advanced LIGO interferometers on 
line in 2014 will extend the distance to 200 Mpc. Short GRBs, if produced by 
mergers as is thought to be the case, predict a cosmological rate density of > 10 
Gpc~^ yr~^ with likely rate of ~300 Gpc~^ yr~^ (O'Shaughnessy et al. 



2005 



O'Shaughnessy et al.||2007l [O'Shaughnessy et al.||2008bl [O'Shaughnessy et al. 



2008a Nakar 2007b). This translates into an Advanced LIGO detection rate 



~10 yr~ . The density estimates include mergers with 7— ray jets not aimed in 
our direction, so not all LIGO detection would be in coincidence with GRBs. 
LIGO is already providing useful upper limits, as with GRB 070201 described in 
Section III. 



6.3 Cosmology 



One of the frontiers of modern cosmology lies at high redshift, -z ^ 6, when the 
first non-linearities developed into gravitationally-bound systems, whose internal 
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evolution gives rise to stars, galaxies, and quasars; and when the light emitted 
from these first collapsed structures diffuse outwards to reionize the Universe. As 



the (temporarily) brightest source of photons in the cosmos (see Figure 15), the 



demise of these first generations of massive stars in GRB explosions defines the 



challenge of elucidating the end of the cosmic "dark ages" (Lamb and Reichart 
2000} [Bromm and Loeb]|2002| ) . 



Apart from revealing a site of high-redshift star-formation, each such high- 
redshift burst has the potential to help constrain local element abundances in 
its host galaxy, information that will be impossible to gather by other means 
until the advent of D ^ 20 m telescopes. Even more exciting, each burst has 
the potential to reveal the extent of intergalactic reionization, at its redshift and 



along that line of sight ( Miralda-Escude 1998 McQuinn et al. 2008 Mesinger 



and Furlanetto 2008). 



At the highest redshifts, z ^ 10, there is growing theoretical evidence indicating 
that the first luminous objects to form were very massive stars M > lOOM©. De- 
pending on whether these stars retain their high masses until death, and whether 
a fast-rotating core is prerequisite to the GRB phenomenon, these Pop III stars 
might provide the progenitors for the most luminous, highest-redshift GRBs 



(Heger et al. 2003). 



The first GRBs and supernovae may also be important for another reason: they 
may generate the first cosmic magnetic fields. Mass loss (e.g. via winds) would 
disperse magnetic flux along with the heavy elements. The ubiquity of heavy 
elements in the Lyman alpha forest indicates that there has been widespread 
diffusion from the sites of these early supernovae, and the magnetic flux could 
have diffused in the same way. This flux, stretched and sheared by bulk motions, 
could be the "seed" for the later amplification processes that generate the larger- 
scale fields pervading disc galaxies. 



6.4 Theoretical Prospects 

What can we expect in the wave of matching theoretical progress? This is more 
difficult to disucss because theory often develops on a shorter timescale than ob- 
servations and experiments, and so we cannot forsee most future developments. 
Although some of the features now observed in GRB sources (especially after- 
glows) were anticipated by theoretical discussions, the recent burst of observa- 
tional discovery has left theory lagging behind. There are, however, some topics 
on which we do believe that there will be steady work of direct relevance to 
interpreting observations. 

One of the most important is the development and use of hydrodynamical codes 
for numerical simulation of GRB sources with detailed physics input. Existing 
two and three dimensional codes have already uncovered some gas-dynamical 



properties of relativistic flows unanticipated by analytical models ( McKinney and 



Blandford 2008 ) but there are some key questions that they cannot yet address. 



In particular, higher resolution is needed because even a tiny mass fraction of 
baryons loading down the outflow severely limits the maximum attainable Lorentz 
factor. What is more, jets are undoutedly susceptible to hydrostatic and hydro- 
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Figure 15 A 360° vista showing the entire sky, with visible structures stretching back in 
distance, time and redshift. The most distant light we observe comes from the radiation 
left over from the Big Bang: the CMB. As we descend the chart, we find the most 
distant objects known, followed by a web of Sloan Digital Sky Survey (SDSS) quasars 
and galaxies. Closer to home, we start to see a collection of familiar "near" galaxies 
(triangles). Also marked are all Swift GRBs with known distances (blue stars); SN 
1997fF, the most distant type la supernova aX z = 1.7; and the archetypal large galaxy 
cluster, the Coma cluster. The redshift distances of most distant GRBs are comparable 
to the most distant galaxies and quasars (adapted from Ramirez- Ruiz ||2006a ) 
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magnetic turbulence. We must wait for useful and affordable three dimensional 
simulations before we can understand the nonlinear development of instabilities. 
Well-resolved three dimensional simulations are becoming increasingly common 
and they rarely fail to surprise us. The symmetry-breaking involved in tran- 
sitioning from two to three dimensions is crucial and can lead to qualitatively 
new phenomena. A particularly important aspect of this would be to link in a 
self-consistent manner the flow within the accretion disk to that in the outflow- 
ing gas, allowing for feedback between the two components. A self-consistent 
model incorporating inflow and outflow must also explain how some fraction of 
the material can acquire more than its share of energy (i.e., a high enthalpy or 
p/p)- 

Particle acceleration and cooling is another problem that seems ripe for a more 
sophisticated treatment. Few would dispute the statement that the photons which 
bring us all our information about the nature of GRBs are the result of particle 



acceleration in relativistic shocks (. 


Uedvedev and Spitkovsky 2008 


Spitkovsky 


2008 


Ramirez-Ruiz et al. 


2007 


Katz et al.j|2007j) or turbulence ( 


Goodman and 


MacFadyen 


2007 


Zhang et al. 


2009 


Couch et al. 20081. Since charged particles 



radiate only when accelerated, one must attempt to deduce from the spectrum 
how the particles are being accelerated, why they are being accelerated, and to 
identify the macroscopic source driving the microphysical acceleration process. 
CoUisionless shocks are among the main agents for accelerating ions as well as 



electrons to high energies whenever sufficient time is available (e.g., Blandford 



and Eichler||1987 Achterberg et al^|2001 1. Particles reflected from the shock and 



from scattering centres behind it in the turbulent compressed region have a good 
chance of experiencing multiple scattering and acceleration by first-order Fermi 
acceleration when coming back across the shock into the turbulent upstream re- 
gion. Second-order or stochastic Fermi acceleration in the broadband turbulence 
downstream of coUisionless shocks will also contribute to acceleration. In addi- 
tion, ions may be trapped at perpendicular shocks. The trapping is a consequence 
of the shock and the Lorentz force exerted on the particle by the magnetic and 
electric fields in the upstream region. With each refiection at the shock the par- 
ticles gyrate parallel to the motional electric field, picking up energy and surfing 
along the shock surface. All these mechanisms are still under investigation, but 
there is evidence that shocks play a most important role in the acceleration of 
cosmic rays and other particles to very high energies. 

Another topic on which further work seems practicable concerns the kinematics 
of ultrarelativistic jets ( Granot||2007 |. While it seems probable that we are using 
the correct ingredients of special relativity and a coUimated outfiow, it is equally 
true that no detailed model yet commands majority of support. We can still 
expect some surprises from studies related to the appearence of relativistic shocks 
in unsteady jets. 

The most interesting problem remains, however, in the nature of the central 
engine and the means of extracting power in a useful coUimated form. In all 
observed cases of relativistic jets, the central object is compact, either a neutron 
star or black hole, and is accreting matter and angular momentum. In addition, in 
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most systems there is direct or indirect evidence that magnetic fields are present 
- detected in the synchrotron radiation in galactic and extragalactic radio sources 
or inferred in collapsing supernova cores from the association of remnants with 
radio pulsars. This combination of magnetic field and rotation may be very 
relevant to the production of relativistic jets. 



7 CONCLUSIONS 



Thanks primarily to the burst discoveries and observations of the Swift satellite, 
the last five years have been tremendously productive ones for GRB research. 
The identification of short burst afterglows has confirmed long-held suspicions 
that GRBs have at least two fundamentally different types of progenitors; subse- 
quent studies of short burst afterglows and host galaxies have furnished valuable 
information on the nature of their progenitors and provided hopeful indications 
for next-generation gravity- wave observatories. The discovery of the first two 
bursts at z > 6, before sources reionized most of the hydrogen in the universe, 
has proven the value of GRBs as probes of the earliest cosmic epochs and and ex- 
tended GRB observations beyond the redshifts of the most distant known quasars 
(see Figure 16). A flood of prompt burst alerts has fed the queues of more than 
a dozen dedicated robotic telescopes, and prompted fast-response multi-epoch 
high-resolution spectroscopy from the largest telescopes. Overall, bright GRBs 
continue to attract the attention of astronomers of all types, with premier facilities 
across the electromagnetic spectrum poised to respond to the next spectacular 
event. In the near future, we hope that gravity- wave and high-energy neutrino 
astronomers will be rewarded for their decades of persistent effort. 

The complexity - not to mention sheer volume - of data in the Swift era have 
inevitably raised challenges to the prior interpretation of bursts and their after- 
glows. However, careful consideration of the biases inherent to the Swift mission 
observing strategies, along with the large number of events now available for 
analysis, are gradually enabling the construction of a single coherent picture, via 
a multidisciplinary approach that addresses data from across the electromagnetic 
spectrum. It is one of the challenges of contemporary research to infer the un- 
derlying physical structure of GRBs in the belief that this is simple, despite the 
complex character of the observations. 

GRBs provide us with an exciting opportunity to study new regimes of physics. 
As we have described, our rationalization of the principal physical considerations 
combines some generally accepted features with some more speculative and con- 
troversial ingredients. When confronted with observations, it seems to accom- 
modate their gross features but fails to provide us with a fully predictive theory. 
What is more valuable, though considerably harder to achieve, is to refine models 
like the ones advocated here to the point of making quantitative predictions, and 
to assemble, assess and interpret observations so as to constrain and refute these 
theories. What we can hope of our present understanding is that it will assist us 
in this endeavour. 

There are good prospects for a continued high rate of discovery going forward. 
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Figure 16 High-z record holders. The history of most distant objects detected in cate- 
gories of galaxies, quasars, and GRBs (R. McMahon, N. Tanvir, priv. comm.). 



Swift is likely to remain operational for many years and Fermi and A GILE, now 
on orbit, are providing new insights at higher energies. Instruments currently 
in development or planning have the potential to enable further qualitative ad- 
vances. GRBs are among the most extraordinary of astronomical phenomena 
and will, with our present and future capabilities, continue providing a unique 
window into the extreme reaches of the Universe. 
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